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ABSTRACT 

We report observations of the linear polarisation of a sample of 50 nearby southern bright stars 
measured to a median sensitivity of ~4.4 xlO^®. We find larger polarisations and more highly 
polarised stars than in the previous PlanetPol survey of northern bright stars. This is attributed 
to a dustier interstellar medium in the mid-plane of the Galaxy, together with a population 
containing more B-type stars leading to more intrinsically polarised stars, as well as using a 
wavelength more sensitive to intrinsic polarisation in late-type giants. Significant polarisation 
had been identified for only six stars in the survey group previously, whereas we are now able 
to deduce intrinsic polarigenic mechanisms for more than twenty. 

The four most highly polarised stars in the sample are the four classical Be stars (a 
Eri, a Col, tj Cen and a Ara). For the three of these objects resolved by interferometry, the 
position angles are consistent with the orientation of the circumstellar disc determined. We 
find significant intrinsic polarisation in most B stars in the sample; amongst these are a number 
of close binaries and an unusual binary debris disk system. However these circumstances do 
not account for the high polarisations of all the B stars in the sample and other polarigenic 
mechanisms are explored. Intrinsic polarisation is also apparent in several late type giants 
which can be attributed to either close, hot circumstellar dust or bright spots in the photosphere 
of these stars. Aside from a handful of notable debris disk systems, the majority of A to K 
type stars show polarisation levels consistent with interstellar polarisation. 

Key words: polarisation - techniques: polarimetric - ISM: magnetic fields - stars: binary - 
stars: emission-line. Be - stars: late-type - stars: giants - debris disks - stars: B spectral type 


1 INTRODUCTION 

The measured linear polarisation of starlight falls into two cate¬ 
gories: it is either intrinsic to the star and its immediate environ¬ 
ment (i.e. the star system), or it is interstellar in origin. Interstellar 
polarisation is the result of interstellar dust particles aligning with 
the Galactic magnetic field. Studies of this polarisation reveal de¬ 
tails of the dust distribution and the Galactic magnetic field struc¬ 
ture (Heiles 1996) as well as the size and nature of the dust particles 
(Whittet et al. 1992; Kim et al. 1994). Linear polarimetry is there¬ 
fore a powerful technique for investigating the interstellar medium. 

Linear polarisation measurements of stars have, until fairly 
recently, been limited to sensitivities in fractional polarisation of 
~10“‘'. At this level most nearby (d < 100 pc) stars are found to be 
unpolarised (Tinbergen 1982; Leroy 1993a,b, 1999). A new genera¬ 
tion of high precision polarimeters has recently been developed that 
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allow polarisation measurements of starlight at the parts per million 
level (Hough et al. 2006; Wiktorowicz & Matthews 2008; Bailey 
et al. 2015; Wiktorowicz & Noli 2015). With such instruments it is 
possible to detect significant interstellar polarisation even in nearby 
stars. Bailey et al. (2010) reported such a survey of 49 bright nearby 
northern hemisphere stars (the PlanetPol survey) and found that the 
majority of stars show significant polarisation, in most cases con¬ 
sistent with an interstellar origin. 

In this paper we present the first survey of the polarisation of 
southern bright nearby stars measured at the parts per million level. 
The survey utilised the newly commissioned HIPPI (High Preci¬ 
sion Polarimetric Instrument) (Bailey et al. 2015) to observe 50 of 
the brightest southern stars within 100 pc. The observing methods 
are described in detail in Section 2. In contrast to the PlanetPol 
survey (Bailey et al. 2010) described above, the southern sample 
reported here contains many stars that are intrinsically polarised. 

Intrinsic polarisation is well known to be present in Be stars as 
a result of scattering from a transient gaseous circumstellar disc and 
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our sample includes four classical Be stars. These are indeed the 
four highest polarisations measured. However, compared to other 
stars we observed we also find higher than average polarisations 
in several B stars, not known to be Be stars. These results are dis¬ 
cussed in Sections 4.8 and 4.9 respectively. Other sources of in¬ 
trinsic polarisation found in the sample are the presence of a de¬ 
bris disk' around the star (Section 4.4), and polarisation, probably 
mostly due to dust scattering, in late type giants (Section 4.10). 

Another area of interest is the level to which normal stars are 
intrinsically polarised. The quiet Sun has been directly measured, 
and shows polarisation of <3 xl0“^ (Kemp et al. 1987). Starspots in 
active stars will show higher polarisation (Strassmeier 2009; Wik- 
torowicz 2009). As well as being interesting in their own right, 
starspots can complicate the growing and exciting field of exoplanet 
polarimetry (Wiktorowicz 2009). 

Polarimetry can be used to detect unresolved hot-Jupiter type 
planets (Seager et al. 2000; Lucas et al. 2006, 2009), as a differ¬ 
ential technique to detect planets in imaging (Schmid et al. 2005; 
Keller 2006) or transit observations (Kostogryz et al. 2015), and to 
characterise exoplanet atmospheres (Bailey 2007; Kopparla et al. 
2014). In each of these applications the signal from the planet is 
small. The effectiveness of the technique will depend upon an un¬ 
derstanding of non-planetary sources of polarisation as much as 
planetary ones. For this reason an investigation and characterisation 
of the primary sources of stellar linear polarisation is important. 


2 OBSERVATIONS 
2.1 The sample stars 

Stars selected for observation were south of the equator, were at a 
distance of less than 100 pc’, and have V magnitude brighter than 
3.0^’. 54 stars met these criteria and we have observed 50 of them. 
The stars observed are listed in Table 2.1. This table gives the V 
magnitude and spectral type as listed in the SIMBAD database"', 
the distance derived from the Hipparcos catalogue parallax (Perry¬ 
man et al. 1997; van Leeuwen 2007), the Galactic coordinates (also 
from the Hipparcos catalogue), and rotational velocity (V sin i) 
from a variety of sources as indicated in Table 2.1. Where a dis¬ 
tinction has been made in the reference, the value given is for the 
primary in multiple systems. The number of companions that fall 
within the aperture was determined in the first instance from Bum- 
ham’s Celestial Handbook (Burnham 1978a,b,c) and from Eggleton 
& Tokovinin (2008), then followed up with a variety of sources as 
indicated in the footnote of Table 2.1. 

Previous polarisation measurements for the stars are listed in 
Table 2 and have been taken from the agglomerated polarisation 
catalogue of Heiles (2000). Only the degree of polarisation is listed. 

* Throughout this paper disk is used in reference to debris disks as pre¬ 
ferred by the debris disk research community, disc being used otherwise. 

^ One star, HIP 89931, has a distance greater than 100 pc as a result of 
revision of its parallax after the sample was selected. 

^ A V magnitude of 3.0 corresponds to an absolute magnitude of -2.0 at 
100 pc. 

HIP 60718, has no listed luminosity class in SIMBAD, and so instead we 
use the spectral type given by Bartkevicius & Gudas (2001). There are a 
number of multiple systems listed, HIP 93506 and HIP 65474, for example, 
are both noteworthy double star systems where both stars are of B spectral 
type, where the combined spectral type is given by SIMBAD this is what is 
given in Table 2.1. If a companion also falls within the HIPPI aperture this 
is marked with an in the succeeding column. 


The position angle can be found in the original catalogue, but for 
almost all the measurements the polarisation is not significant and 
the position angle is therefore meaningless. The only stars with sig¬ 
nificant polarisations in this catalogue are the four Be stars, which 
will be discussed later, the SRB (semi-regular) type pulsating vari¬ 
able HIP 112122 ip Gru) which has a listed polarisation angle of 
155.2°±4.6, and possibly the short period binary and variable of 
the p Cep type HIP 65474 (Spica) with a listed polarisation angle 
of 147.0°±0.0. 

The Heiles (2000) catalogue does not include the slightly more 
accurate previous polarisation measurements of nearby stars made 
by Tinbergen (1982). This study includes a number of the stars in 
our sample and is therefore listed separately in Table 2. The mea¬ 
surements of Tinbergen were made in three colour bands. We have 
tabulated the averaged band I and II measurements. Of the stars ob¬ 
served by Tinbergen (1982) only HIP 89931 (S Sgr) and Spica have 
a significant polarisation. 

In addition to the measurement listed in Table 2, it might be 
noted that Sirius (HIP 32349) also has a reported measurement in 
Bailey et al. (2015) where it was used as a low polarisation standard 
to determine telescope polarisation. When the telescope polarisa¬ 
tion is subtracted the polarisation of Sirius is 3±4 ppm. 

2.2 Observation methods 

The observations were obtained with the HIPPI (High Precision 
Polarimetric Instrument) polarimeter (Bailey et al. 2015) over the 
course of four observing runs (May 8th to May 12th 2014, August 
28th to September 2nd 2014, May 21st to 26th 2015 and June 25th 
to 29th 2015) at the 3.9 m Anglo Australian Telescope (AAT). The 
AAT is located at Siding Spring Observatory near Coonabarabran 
in New South Wales, Australia. HIPPI was mounted at the f/8 
Cassegrain focus of the telescope where it had an aperture size of 
6.7". The seeing varied between ~ 1.5" and > 4" with a typical see¬ 
ing around 2-2.5". 

HIPPI is a high precision polarimeter, with a sensitivity in 
fractional polarisation of ~4.3 xl0“® on stars of low polarisation 
and a precision of better than 0.01% on highly polarised stars (Bai¬ 
ley et al. 2015). It achieves this by the use of a Ferroelectric Liq¬ 
uid Crystal (FLC) modulator operating at a frequency of 500 Hz 
to eliminate the effects of variability in the atmosphere. From the 
modulator the light passes through a Wollaston prism that acts as 
a beam splitter, directing the light into two Photo Multiplier Tube 
(PMT) detectors. Second stage chopping, to reduce systematic ef¬ 
fects, is accomplished by rotating the entire back half of the in¬ 
strument after the filter wheel through 90° in an ABBA pattern, 
with a frequency of once per 40 seconds. An observation of this 
type measures only one Stokes parameter of linear polarisation. To 
get the orthogonal Stokes parameter the entire instrument is rotated 
through 45° and the sequence repeated. The rotation is performed 
using the AAT’s Cassegrain instrument rotator. In practice we also 
repeat the observations at telescope position angles of 90° and 135° 
to allow removal of instrumental polarisation. The stars in the sur¬ 
vey were observed for a total integration time of 12 minutes in each 
(Q and U) Stokes parameter per observation. This gave a median 
sensitivity in fractional polarisation of ~4.4 xl0“®. The observing, 
calibration and data reduction methods are described in full detail 
in Bailey et al. (2015). 

Typically, a sky measurement was acquired at each telescope 
position angle an object was observed in, and subtracted from the 
measurement. The duration of the sky measurements was 3 min¬ 
utes per Stokes parameter. For particularly bright objects (or those 
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Table 1. Properties of survey stars. 


HIP 

BS 

Other Names 

V 

Spectral 

Comp 

Dist 

RA 

Dec 

Galactic 

V sin i 

f-IR/Ft 




mag 

Type 

in Ap“ 

(pc) 

(hh mm) 

(dd mm) 

Long 

Eat 

(km/s) 

(xlO-®) 

2021 

98 

ySHyi 

2.79 

GOV 


7.5 

00 26 

-77 15 

304.77 

-39.78 

4.1 ± 

IQdM 


2081 

99 

a Phe 

2.37 

K0.5IIIb 


25.9 

00 26 

-42 18 

320.00 

-73.98 

< 



3419 

188 

P Get 

2.01 

KOIII 


29.5 

00 44 

-17 59 

111.33 

-80.68 

3.6 ± 



7588 

472 

a Eri, Achemar 

0.46 

B6Vep 


42.9 

01 38 

-57 14 

290.84 

-58.79 

207 ± 

9’' 

0.027 

9236 

591 

a Hyi 

2.84 

FOIV 


22.0 

01 59 

-61 34 

289.44 

-53.76 

118 ± 

12Z« 

34.24 

13847 

897 

eEri 

2.90 

A4III 


49.4 

02 58 

-40 18 

247.86 

-60.74 

70 ± 17^ 


18543 

1231 

y Eri 

2.94 

Ml mb 


62.3 

03 58 

-13 30 

205.16 

-44.47 

3.8 ± 

0.3'’^ 


26634 

1956 

a Col 

2.65 

B9Ve 


80.1 

05 40 

-34 04 

238.81 

-28.86 

184 ± 

5’' 

1.484 

30438 

2326 

a Car, Canopus 

-0.74 

A9II 


94.8 

06 24 

-52 42 

261.21 

-25.29 

9.0 ± 

1.0® 

0.033 

32349 

2491 

a CMa, Sirius 

-1.46 

AIV 


2.6 

06 45 

-16 43 

227.23 

-08.89 

16 ± 

1® 


39757 

3185 

p Pup 

2.81 

E5II 


19.5 

08 08 

-14 18 

243.15 

04.40 

15 ± 

jGG 

160.7 

42913 

3485 

(5Vel 

1.95 

AlVa(n) 

* 

24.7 

08 45 

-54 43 

272.08 

-07.37 

150 ± 15^^ 

4.567 

45238 

3685 

P Car 

1.69 

Aim 


34.7 

09 13 

-69 43 

285.98 

-14.41 

146 ± 

2®’ 

6.700 

46390 

3748 

a Hya, Alphard 

1.97 

K3ii-m 


54.0 

09 28 

-08 40 

241.49 

29.05 

< 

1.0®“ 

3.359 

52727 

4216 

p Vel 

2.69 

G6m 

* 

35.9 

10 47 

-49 25 

283.03 

08.57 

6.4 ± 

qJMM 

19.94 

59803 

4662 

y Crv 

2.58 

B8m 


47.1 

12 16 

-17 33 

290.98 

44.50 

30 ± 

9”' 


60718 


a Cm 

0.81 

B0.5IV 


98.7 

12 27 

-63 06 

300.13 

-00.36 

88 ± 

5GG 

0.001 

60965 

4757 

(5 Crv 

2.94 

AOIV(n) 


26.6 

12 30 

-1631 

295.47 

46.04 

236 ± 

24Z® 


61084 

4763 

y Cm 

1.64 

M3.5m 


27.2 

1231 

-57 07 

300.17 

05.65 

4.7 ± 

0.9® 


61359 

4786 

P Crv 

2.64 

G5II 


44.7 

12 34 

-23 24 

297.87 

39.31 

4.2 ± 

j QdMM 

5.154 

61585 

4798 

a Mus 

2.68 

B2IV-V 


96.7 

12 37 

-69 08 

301.66 

-06.30 

114 ±11*^ 

0.179 

61932 

4819 

y Cen 

2.17 

AlIV+AOIV 

* 

39.9 

12 42 

-48 58 

301.25 

13.88 

79 ± 

3GG 


61941 


y Vir 

2.74 

EOIV+EOIV 

* 

11.7 

12 42 

-01 27 

297.83 

61.33 

36 ± 

4ZR 


64962 

5020 

y Hya 

3.00 

G8m 


41.0 

13 19 

-23 10 

311.10 

39.26 

3.7 ± 

I.0®“ 

2.931 

65109 

5028 

L Cen 

2.73 

A2V 


18.0 

13 20 

-36 43 

309.41 

25.79 

90 ± 

2D 

14 

65474 

5056 

a Vir, Spica 

0.97 

BlIII-IV 

* 

77.0 

13 25 

-11 10 

316.11 

50.84 

140 ± 

6GG 

0.003 

68933 

5288 

6 Cen 

2.05 

KOIII 


18.0 

14 07 

-36 22 

319.45 

24.08 

1 2CC 



71352 

5440 

r} Cen 

2.31 

B1.5Vne 


93.7 

14 36 

-42 09 

322.77 

16.67 

291 ± 32’' 

0.443 

71683 

5459 

a Cen 

0.10 

G2V+K1V 

* 

1.3 

14 40 

-60 50 

315.73 

-00.68 

2.3 ± 

0 . 5 “'^ 


72622 

5531 

a~ Lib 

2.75 

A3V 

* 

23.2 

14 50 

-16 02 

340.32 

38.01 

102 ± 

IqZ® 

16.25 

74785 

5685 

ySLib 

2.62 

B8Vn 

* 

56.8 

15 17 

-09 22 

352.02 

39.23 

250 ± 

24 ^ 

3.652 

74946 

5671 

y TrA 

2.89 

AIV 


56.4 

15 19 

-68 41 

315.71 

-09.55 

199 ± 20 ^'^ 

5.391 

77952 

5897 

ySTrA 

2.85 

ElV 


12.4 

15 55 

-63 26 

321.85 

-07.52 

^ 2 ^Ma 



79593 

6056 

d Oph 

2.75 

M0.5m 


52.5 

15 14 

-03 41 

8.84 

32.20 

7.0 ± 

0.0“ 


82396 

6241 

6 ScO 

2.29 

Kim 


19.5 

16 50 

-34 18 

348.81 

06.56 

< 

1.0®“ 


84012 

6378 

rj Oph 

2.42 

A2IV-V 

* 

27.1 

17 10 

-15 44 

6.72 

14.01 

23 ± 

2Z® 

19.01 

85792 

6510 

a Ara 

2.95 

B2Vne 


82.0 

17 32 

-49 52 

340.75 

-08.83 

269 ± 

13’' 

15.62 

86228 

6553 

6 Sco 

1.86 

El III 


92.1 

17 37 

-43 00 

347.14 

-05.98 

125’'® 


3.023 

88635 

6746 

y- Sgr 

2.99 

Kim 


29.7 

18 06 

-30 25 

0.92 

-04.54 

< 

1.0®“ 

146.7 

89931 

6859 

<5 Sgr 

2.71 

K3ma 


106.6 

1821 

-29 50 

3.00 

-70.15 

3.6 ± 

1.0®“ 

1.305 

90185 

6879 

eSgr 

1.85 

B9.5m 

* 

43.9 

18 24 

-34 23 

359.19 

-09.81 

175”® 


31 

90496 

6913 

A Sgr 

2.81 

KOIV 


24.0 

18 28 

-25 25 

7.66 

-06.52 

< 

2 ^dM 


92855 

7121 

o-Sgr 

2.06 

B2V 


69.8 

18 55 

-26 18 

9.56 

-12.43 

205”® 


0.025 

93506 

7194 

^Sgr 

2.61 

A2.5Va 

* 

27.0 

19 03 

-29 53 

6.84 

-15.35 

77 ± 

8Z® 

7.111 

100751 

7790 

a Pav 

1.91 

B2IV 

* 

54.8 

20 26 

-56 44 

340.9 

-35.19 

15 ± 

3GG 

0.038 

107556 

8322 

d Cap 

2.83 

A7m 

* 

11.9 

21 47 

-16 08 

37.60 

-46.01 

93 ± 

lOGG 


109268 

8425 

a Gru 

1.71 

B6V 


31.0 

22 08 

-46 58 

349.99 

-52.47 

259 ± 26^^® 


110130 

8502 

a Tuc 

2.82 

K3m 

* 

61.3 

22 19 

-60 16 

330.22 

-47.96 

1.9 ± 

1.3®“ 

6.579 

112122 

8636 

P Gru 

2.11 

M5m 


54.3 

22 43 

-46 53 

346.27 

-57.95 

2.5 ± 

1.5’’ 


113368 

8728 

a PsA, Fomalhaut 

1.16 

A4V 


7.7 

22 58 

-29 37 

20.49 

-64.91 

93 ± 

^ZR 

22 


a - An asterisks indicates a companion also lies within the aperture. Sources for these determinations are: Burnham (1978a,b,c), Eggleton & Tokovinin 
(2008), Pourbaix et al. (2004); Spica - Harrington et al. (2009); e Sgr - Hubrig et al. (2001); ji Lib - Roberts et al. (2007). 

b - Fractional infrared excess from McDonald et al. (2012), except for l Cen and Fomalhaut (Marshall et al. 2015), and e Sgr for which we derived the 
fractional luminosity based on excesses reported by Rhee et al. (2007) and Chen et al. (2014). 
c - There are at least two stars within the aperture, there maybe as many as four (see Section 3.1.2). 

Rotational velocity references: dM - De Medeiros et al. (2014), Y - Yudin (2001), ZR - Zorec & Royer (2012), R - Royer et al. (2002), M - Massarotti et al. 
(2008), dMM - de Medeiros & Mayor (1999), A - Abt et al. (2002), GG - Gl^bocki & Gnaciriski (2005), C - Cummings (1998), W - Wollf et al. (2007), D - 
Diaz et al. (2011), We - Weise et al. (2010), Ma - Mallik et al. (2003), vB - van Belle (2012), Z - Zamanov et al. (2008). 
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Table 2. Previous polarisation measurements. 


HIP 

Polarisation (Heiles) 
p(%) 

Polarisation (Tinbergen)" 
q (ppm) u (ppm) 

2021 

0.050±0.035 

60 ±60 

0 ± 60 

2081 

0.011±0.006 

60 ±60 

-70 ± 60 

3419 

0.004±0.005 

40 ±60 

-60 ± 60 

7588'’ 

0 .011±0.002 

10 ±60 

-10± 60 

9236 


30 ±60 

-10± 60 

13847 

0.003±0.003 



26634'’ 

0.150±0.100 



32439 


0±60 

-40 ± 60 

46390 

0.150±0.120 

30 ±60 

10 ± 60 

52727 


30 ±60 

-40 ± 60 

59803 

0.050±0.120 



60965 

0.030±0.120 



61359 


120 ± 60 

-30 ± 60 

61941 

0.030±0.120 

60 ±60 

70± 60 

65109 


-30 ± 60 

110 ± 60 

65474^' 

0 .010±0.000 

-220 ± 60 

-10± 60 



-150 ±60 

-60 ± 60 

71352* 

0.040±0.100 



71683 


-70 ± 60 

30± 60 

74785 

0.050±0.120 



77952 


-40 ± 60 

40± 60 

79593 

0 .020±0.120 



82396 


60 ±60 

-90 ± 60 

84012 

0.026±0.004 



85792* 

0.610±0.035 



86228 

0.030±0.035 



89931 


-110 ±60 

-520 ± 60 

90496 


90 ±60 

-150± 60 

107556 


-40 ± 60 

-40 ± 60 

109268 

0.003±0.005 

-70 ± 60 

10 ± 60 

110130 

0.011±0.004 



112122 

0.031±0.005 



113368 

0.005±0.005 

10 ±60 

-110± 60 


a - The polarisations were originally reported with units of 10“^ 
b - Classical Be star. 

c - Spica appears in Tinbergen’s table twice; observed from 
Hartebeesportdam (top) and LaSilla (bottom). 


known to be highly polarised) observed under a moonless sky, the 
sky signal is negligible and a dark measurement was sufficient for 
calibration purposes. Dark measurements are obtained with the use 
of a blank in the filter wheel and the instrument otherwise config¬ 
ured identically. Whether a sky subtraction or only a dark subtrac¬ 
tion was carried out is indicated in Table 5. These subtractions were 
carried out as the first step of the data reduction routine. 

During the observations at the AAT, zero point calibration was 
carried out by reference to the average of a set of observed stars that 
have been measured as having negligible polarisation. This method 
was used owing to a lack of other available unpolarised standard 
stars in the southern hemisphere, and allowed us to determine the 
telescope polarisation for the equatorially mounted AAT. Using this 
method the polarisation of the telescope was determined as 48±5 
xl0“® during the Aug-Sep run (Bailey et al. 2015). This value was 
used for correction of both 2014 runs, with the uncertainty in the 
measurement contributing to the stated errors. Shortly before the 
May 2015 run the telescope mirror was re-aluminised necessitating 
a new determination of the telescope polarisation. This procedure 
was carried out as before, with the details given in Table 3. HIPPI 


Table 3. Low polarisation star measurements to determine telescope polar¬ 
isation (TP) for the May and June 2015 runs in the g' filter. The adopted TP 
for June 2015 uses all seven measurements, that for May 2015 just those 
acquired in that month. 


Star 

Date 

P (ppm) 

en 

Sirius 

23 May 

39.2 ± 1.4 

91.0 ± 2.1 


24 May 

40.2 ± 0.8 

89.5 ± 1.1 

BS 5854 

22 May 

30.4 ± 5.1 

70.2 ± 9.7 


26 May 

37.6 ± 4.0 

84.4 ± 6.1 

Adopted TP 

May 2015 

35.5 ± 1.4 

84.7 ± 2.3 

BS 5854 

26 Jun 

43.3 ± 4.4 

93.4 ± 5.9 


27 Jun 

32.9 ± 3.9 

70.6 ± 6.8 

Leo 

26 Jun 

43.6 ± 2.4 

87.0 ± 3.2 

Adopted TP 

Jun 2015 

36.5 ± 1.2 

84.8 ± 1.4 


was left in place at the AAT Cassegrain focus between the May and 
June 2015 runs. Three additional measurements of low polarisation 
standards (fi Leo, and two of BS 5854) were made in June. To in¬ 
crease the precision of the determination for June 2015 the average 
of all the May and June measurements was adopted as the telescope 
polarisation for that run - this is reasonable given the negligible dif¬ 
ference in the determinations. The telescope polarisation was found 
to be 36.5±1.2xl0-^ 

A number of stars with known high polarisations (~l-5%) 
were observed during the Aug-Sep run, and used to help calibrate 
the modulation efficiency of the telescope as reported by Bailey 
et al. (2015). The same stars were used to determine the position 
angle zero-point for the Aug-Sep run, whereas a different set were 
needed for the other runs; these were HD 147084, HD 154445, HD 
160529 and HD 187929 in May 2014; HD 147084 and HD 154445 
in May 2015; and HD 147084 and HD 154445 in June 2015. The 
precision of each determination is less than 1 degree, based on the 
consistency of the calibration provided by the different reference 
stars which themselves have uncertainties of this order. 

HIPPI has been tuned to be most sensitive at the blue end of 
the visible spectrum, in order to investigate polarisation induced in 
exoplanet atmospheres by Rayleigh scattering from clouds (Bailey 
et al. 2015). This it achieves through the use of PMTs with ultra¬ 
bialkali photocathodes, giving a quantum efficiency of 43% at 400 
nm. HIPPI is equipped with a number of filters. The initial obser¬ 
vations reported here were carried out exclusively using the SDSS 
g' (400-550 nm) filter. 

The g' band is centred on 475 nm and is 150 nm in width, 
which results in the precise effective wavelength^ and modulation 
efficiency changing with star colour. Table 4 gives the effective 
wavelength and modulation efficiency for various spectral types 
based on a bandpass model as described in Bailey et al. (2015). 
It combines the attenuation of the atmosphere, as well as PMT re¬ 
sponse and other instrumental responses as well as that of the g' 
filter. No interstellar extinction (i.e. E(B - V) = 0) has been ap¬ 
plied here. This is reasonable as at 100 pc the extinction is well 
below that which can actually be measured by photometry (Bailey 
et al. 2010). The efficiencies in Table 4 were applied to the raw 
Qjl and (/// polarisations of each object according to its spectral 


^ Effective wavelength is the intensity weighted mean of the wavelength of 
observed photons. 
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Table 4. Effective wavelength and modulation efficiency for different spec¬ 
tral types according to bandpass model. 


Spectral 

Type 

Effective 
wavelength (nm) 

Modulation 
efficiency (%) 

BO 

459.1 

87.7 

AO 

462.2 

88.6 

FO 

466.2 

89.6 

GO 

470.7 

90.6 

KO 

474.4 

91.6 

MO 

477.5 

92.0 

M5 

477.3 

91.7 


type given in Table 2.1 (we apply a linear interpolation between the 
given types). 

Follow up observations were made for one object using SDSS 
r' and 425 nm short pass (425SP) filters on 24/5/2015. The same 
procedures were used for calibrating telescope polarisation, angu¬ 
lar calibration and efficiency as for the observations made using the 
SDSS g' band. The r' filter is 150 nm in width; it is centred on 625 
nm but HIPPFs PMTs decrease in efficiency steadily from ~525 
nm, giving a shorter effective wavelength (this is depicted graphi¬ 
cally in Bailey et al. (2015)). The response of the 425SP filter at the 
blue end is truncated by the detector response at ~350 nm. 

2.3 Accounting for patchy cloud 

Some of the observations were made in cloudy conditions. For three 
objects (HIP 88635, HIP 107556 and to a lesser degree HIP 90496) 
the precision is considerably worse than the median (see Table 5) 
owing to particularly poor weather on September 1st, 2014. Thick, 
rapidly moving patchy (mostly nimbus) cloud was constant for 
much of the night, and seeing was worse than 4" at one point. Sim¬ 
ilarly patchy cloud also affected observations on September 2nd, 
but not to the same degree. 

To account for this we removed the most cloud affected parts 
of the observations. The HIPPI observing technique allows for the 
efficient flagging and elimination of bad data. The data are taken 
in lots of 20 one second integrations per rotation. The Stokes pa¬ 
rameters Q or U and I are determined for each integration using 
a Mueller matrix method (see Bailey et al. (2015) for more de¬ 
tails). The errors are then determined based on the number statis¬ 
tics of all the integrations in a rotation. Integrations with a Stokes 
1 determination of less than 5% of the maximum observed for the 
corresponding Q or U Stokes parameter were discarded as cloud 
affected. The precision of the measurements is improved by remov¬ 
ing the cloud affected integrations. However, each observation was 
initiated in as clear as possible conditions to ensure an appropriate 
maximum 1. 

The 5% threshold was arrived at from very many repeated ob¬ 
servations of HIP 2081 on September 1st, producing hundreds of 
integrations, made during the most varied observing conditions. By 
examining a moving average of one Stokes parameter with /, it was 
found that the accuracy of polarisation determinations was not af¬ 
fected, only the precision, at least down to a threshold of 3%. 


3 RESULTS 

The resulting polarisation measurements for the 50 stars are given 
in Table 5. This table lists the normalized Stokes parameters q = 


Qjl and u = Ujl, on the equatorial system, and the degree of po¬ 
larisation and position angle obtained by combining the q and u 
measurements as p = fq^ + u^. The polarisations and Stokes pa¬ 
rameters are in units of parts per million (ppm, equal to 10“'’) in 
fractional polarisation. 

The errors quoted are derived from the internal statistics of the 
individual data points included in each measurement as described 
in Section 2.3 and by Bailey et al. (2015). 

3.1 Uncertain results for binaries with aperture scale 
separations 

The process of centring a star in HIPPFs 6.7" aperture is carried 
out by manual scanning to maximise the total signal received by 
the instrument. This process is made difficult when stars in a bi¬ 
nary system have a separation similar to the radius of the aperture, 
particularly if they have a similar apparent magnitude, or if the see¬ 
ing is poor. This can result in the system being off-centre in the 
aperture, and a partial contribution from the secondary. In such in¬ 
stances a small instrumental polarisation is induced that would be 
difficult to calibrate for. 

3.1.1 a Cen 

This particular difficulty was apparent when attempting to acquire 
a Cen A (HIP 71683). The separation of a Cen A (G2V) and B 
(KIV) at the time of our observations was around 4" (the separa¬ 
tion is depicted graphically by Burnham (1978a), and the updated 
parameters given by Pourbaix et al. (2002) are little different). We 
certainly have a significant contribution from component B also, 
and so we have reported the polarisation for the a Cen system as a 
whole. In the discussion that follows it will become clear that the 
degree of polarisation of a Cen is anomalously high when one con¬ 
siders its spectral type and proximity to the Sun, and we ascribe 
this to aperture scale separation of a Cen B inducing instrumental 
polarisation. 

3.1.2 a Cm 

Although not apparent at the time, it is possible that a Cru (HIP 
60718) is affected in the same way. The components a* (B0.5IV) 
and (B1V) are also separated by ~4" (Burnham 1978c; Pourbaix 
et al. 2004). a Cru shows higher polarisation than any other (non- 
Be) B type star in our survey which suggests that its measurement 
is spurious. However, there are other factors that may result in a 
high degree of polarisation for a Cru. It is the star system in the 
survey with the earliest spectral type, and the or* component is itself 
a binary with a separation of ~1 AU, while may also be a binary 
(Burnham 1978c; Pourbaix et al. 2004). These factors are discussed 
with reference to other stars in Sections 4.9 and 4.7. 


4 DISCUSSION 

The discussion begins with a comparison to previous results (4.1) 
and a look at the statistics of the survey (4.2). Thereafter it is di¬ 
vided up roughly by stellar spectral type. Polarigenic mechanisms 
for A-K stars begin with the spatial distribution of polarisation due 
to the interstellar medium (4.3), we then look at debris disk sys¬ 
tems (4.4), Ap stars (4.5) and eclipsing binaries (4.6). B spectral 
types are examined next beginning with close binaries (4.7), then 
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Table 5. HIPPI linear polarisation measurements. 


Star 

Datefs)” 

Cal* 

q (ppm) 

u (ppm) 

P (ppm) 

6 0 

HIP 2021 

28-31/8/14 

D 

-8.6 ± 

2.5 

-1.6 ± 

2.5 

8.8 ± 

2.5 

95.1 ± 

16.3 

HIP 2081 

1/9/14 

S 

-5.4 ± 

4.6 

-82.0 ± 

3.9 

82.1 ± 

4.3 

133.1 ± 

3.2 

HIP 3419 

2/9/14 

S 

-23.1 ± 

8.0 

-4.1 ± 

7.8 

23.5 ± 

7.9 

95.1 ± 

19.1 

HIP 7588 

2/9/14, 24/5/15 

D,S 

969.6 ± 

4.0 

1920.9 ± 

3.6 

2151.8 ± 

3.8 

31.6± 

0.1 

HIP 9236 

2/9/14 

D 

31.7± 

8.2 

-28.1 ± 

8.7 

42.4 ± 

8.4 

159.2 ± 

11.5 

HIP 13847 

2/9/14 

D 

33.6 ± 

8.8 

65.9 ± 

8.0 

74.0 ± 

8.4 

31.5 ± 

6.7 

HIP 18543 

31/8/14,2/9/14 

D 

42.0 ± 

6.7 

2.9 ± 

6.0 

42.1 ± 

6.3 

1.9 ± 

8.2 

HIP 26634 

31/8/14 

S 

-668.4 ± 

6.5 

-273.1 ± 

6.3 

722.0 ± 

6.4 

101.1 ± 

0.5 

HIP 30438 

28-30/8/14, 2/9/14 

D 

-68.9 ± 

1.7 

-89.2 ± 

1.6 

112.8 ± 

1.7 

116.2 ± 

0.9 

HIP 32349^ 

31/8/14,2/9/14, 23-24/5/15 

D 

-3.7 ± 

1.7 

-4.0 ± 

1.7 

5.5 ± 

1.7 

113.8± 

17.8 

HIP 39757 

25/5/15 

S 

-13.4 ± 

5.8 

-12.5 ± 

5.8 

18.3 ± 

5.8 

111.5 ± 

18.2 

HIP 42913 

25/5/15 

S 

-10.9 ± 

8.0 

-43.1 ± 

8.0 

44.5 ± 

8.0 

127.9 ± 

10.4 

HIP 45238 

23/5/15 

s 

-20.5 ± 

3.3 

12.4 ± 

3.4 

23.9 ± 

3.4 

74.4 ± 

8.1 

HIP 46390 

24/5/15 

s 

8.4 ± 

4.7 

-2.8 ± 

4.7 

8.8 ± 

4.7 

170.8 ± 

30.3 

HIP 52727 

12/5/14 

s 

-12.7 ± 

5.3 

-30.4 ± 

4.6 

33.0 ± 

4.9 

123.7 ± 

9.0 

HIP 59803 

12/5/14 

s 

9.7 ± 

4.3 

-68.2 ± 

3.4 

68.9 ± 

3.9 

139.1 ± 

3.5 

HIP 60718 

12/5/14, 26/6/15 

s 

-248.1 ± 

2.3 

-258.7 ± 

1.9 

358.4 ± 

2.1 

113. I ± 

0.3 

HIP 60965 

24/5/14 

s 

-14.3 ± 

4.8 

12.1 ± 

4.8 

18.7 ± 

4.8 

70.0 ± 

14.7 

HIP 61084 

12/5/14 

s 

-19.8 ± 

4.2 

-39.8 ± 

3.5 

44.5 ± 

3.9 

121.8 ± 

5.3 

HIP 61359 

12/5/14 

s 

-29.1 ± 

4.8 

15.3 ± 

4.1 

32.9 ± 

4.4 

76.1 ± 

7.4 

HIP 61585 

12/5/14 

s 

145.7 ± 

4.7 

3.3 ± 

3.4 

145.7 ± 

4.0 

0.7 ± 

1.3 

HIP 61932 

24/5/15 

s 

-37.7 ± 

5.3 

48.8 ± 

5.1 

61.6 ± 

5.2 

63.8 ± 

4.9 

HIP 61941 

24/5/15 

s 

-5.3 ± 

5.0 

-5.4 ± 

5.1 

7.6 ± 

5.0 

112.5 ± 

38.2 

HIP 64962 

12/5/14 

s 

2.2 ± 

5.1 

5.6 ± 

4.5 

6.1 ± 

4.8 

34.3 ± 

47.9 

HIP 65109 

2/9/14 

s 

-25.7 ± 

4.9 

22.2 ± 

4.7 

34.0 ± 

4.8 

69.6 ± 

8.0 

HIP 65474 

24/5/15, 2x29/6/15 

s 

-200.2 ± 

2.0 

42.7 ± 

2.0 

204.7 ± 

2.0 

84.0 ± 

0.6 

HIP 68933 

24/5/15 

s 

-16.9 ± 

4.3 

39.1 ± 

4.3 

42.6 ± 

4.3 

56.7 ± 

5.8 

HIP 71352 

12/5/14 

s 

5987.4 ± 

3.6 

-1757.4 ± 

3.1 

6240.0 ± 

3.4 

171.8 ± 

0.0 

HIP 71683 

12/5/14 

s 

30.3 ± 

3.2 

-30.0 ± 

2.1 

42.6 ± 

2.7 

157.6 ± 

3.7 

HIP 72622 

12/5/14 

s 

7.1 ± 

4.2 

-26.6 ± 

3.6 

27.5 ± 

3.9 

142.5 ± 

8.6 

HIP 74785 

12/5/14 

s 

4.9 ± 

4.1 

150.5 ± 

3.4 

150.6 ± 

3.7 

44.1 ± 

1.6 

HIP 74946 

12/5/14 

s 

3.3 ± 

4.8 

-29.4 ± 

4.1 

29.5 ± 

4.4 

138.2 ± 

9.2 

HIP 77952 

9/2/14, 22/5/15 

s 

5.3 ± 

4.4 

3.7 ± 

4.2 

6.4 ± 

4.3 

17.5 ± 

37.9 

HIP 79593 

12/5/14, 24/5/15 

s 

-50.1 ± 

4.7 

-480.1 ± 

4.2 

482.7 ± 

4.5 

132.0 ± 

0.6 

HIP 82396 

1/9/14 

s 

10.5 ± 

6.6 

27.0 ± 

9.2 

28.9 ± 

7.9 

34.4 ± 

13.8 

HIP 84012 

12/5/14 

s 

25.2 ± 

4.4 

-50.7 ± 

3.8 

56.6 ± 

4.1 

148.2 ± 

4.4 

HIP 85792 

12/5/14 

s 

5113.8± 

4.0 

-1469.1 ± 

3.5 

5320.6 ± 

3.8 

172.0 ± 

0.0 

HIP 86228 

12/5/14 

s 

-149.3 ± 

3.6 

-21.3 ± 

3.0 

150.8 ± 

3.3 

94.1 ± 

1.2 

HIP 88635 

1/9/14 

s 

37.6 ± 

20.3 

-8.3 ± 14.7 

38.5 ± 

17.5 

173.8 ± 

22.4 

HIP 89931 

12/5/14 

s 

-313.1 ± 

5.1 

-502.5 ± 

4.5 

592.1 ± 

4.8 

119.0 ± 

0.5 

HIP 90185 

1/9/14 

s 

38.8 ± 

4.5 

158.2 ± 

4.3 

162.9 ± 

4.4 

38.1 ± 

1.6 

HIP 90496 

1/9/14 

s 

9.6 ± 

10.0 

-53.4 ± 

9.1 

54.2 ± 

9.5 

140.1 ± 

10.5 

HIP 92855 

1/9/14, 22/5/15 

s 

-30.7 ± 

3.5 

-166.9 ± 

3.4 

169.7 ± 

3.4 

129.8 ± 

1.2 

HIP 93506 

1/9/14 

s 

0.9 ± 

4.6 

-28.2 ± 

4.7 

28.2 ± 

4.6 

135.9 ± 

9.3 

HIP 100751 

1/9/14 

s 

-5.7 ± 

4.8 

-85.5 ± 

4.1 

85.6 ± 

4.4 

133.1 ± 

3.2 

HIP 107556 

1/9/14 

s 

2.4 ± 

19.8 

-29.5 ± 

15.6 

29.6 ± 17.7 

137.3 ± 

38.3 

HIP 109268 

31/8/14 

D 

-91.4 ± 

3.4 

21.1 ± 

3.1 

93.8 ± 

3.3 

83.5 ± 

1.9 

HIP 110130 

31/8/14, 26/6/15 

D,S 

-107.4 ± 

1.8 

-74.7 ± 

4.2 

131.1 ± 

4.2 

107.4 ± 

1.8 

HIP 112122 

1/9/14 

s 

330.7 ± 

5.8 

560.2 ± 

8.1 

650.5 ± 

7.0 

29.7 ± 

0.6 

HIP 113368 

28/8/14 

D 

-17.8 ± 

3.2 

-16.6 ± 

3.0 

24.3 ± 

3.1 

111.5 ± 

7.4 


a - Hyphenation indicates object was observed once per day inclusive, 
b - Calibration type: full sky subtraction (S), or dark subtraction only (D). 
c - Sirius was used as a low polarisation standard. 

Be stars (4.8) before concluding with a look at the remainder of B 
stars (4.9). Finally we consider the late giants in the survey (4.10). 

4.1 Comparison with previous observations 

Of the stars highlighted in Section 2.1 as having previous signifi¬ 
cant polarisation measurements, all have significant measurements 


in our survey also. Our determinations are all of the same order of 
magnitude as the previous ones. However, all are significantly dif¬ 
ferent. This is reflective either of the improved precision of HIPPI 
over older instruments, or these three stars are variable in polarisa¬ 
tion - an indication that the polarisation is intrinsic. A less likely 
possibility is that the variation represents the movement of the stars 
with respect to the patchy interstellar medium, but the variability 
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Table 6. Polarisation for spectral classes. 


Spectral 

Class 

N 

HIPPr 

Mean 

d(pc) 

Median 

P (ppm) 

N 

PlanetPol* 

Mean Median 
d (pc) p (ppm) 

B'^ 

9 

65.0 

145.7 

5 

49.6 

36.8 

A 

14 

31.7 

29.6 

18 

28.1 

8.4 

F/G 

10 

45.5 

32.9 

9 

23.5 

9.4 

K 

9 

40.9 

54.2 

12 

38.7 

9.8 

M 

4 

49.1 

263.6 

4 

78.6 

109.2 


a - This work, 
b - Bailey et al. (2010). 

c - Be stars not included for HIPPI; if included: 13, 67.3, 166.6. 


we see would correspond to the movement of many tens of parsecs 
worth of dust based on the polarisation with distance relation found 
by Bailey et al. (2010) for northern stars. 

The best match for previous measurements is with that of Tin¬ 
bergen (1982) for Spica, where we are in agreement within the error 
for both previous measurements in q and significantly different to 
less than 2a in u. The measurement tabulated by Heiles (2000) is 
half the degree of polarisation and 63° different. Spica will be dis¬ 
cussed in detail in Section 4.7.1. 

For the late Giant 6 Sgr we are in agreement with Tinbergen 
(1982) in u, but have a thrice greater q measurement that is sig¬ 
nificant to more than 3cr. For f} Gru we have double the polarisa¬ 
tion of Fleiles (2000) and a polarisation angle that is 126° different, 
again this is significant to more than 3a. This is strong evidence for 
variable intrinsic polarisation. These two stars will be discussed in 
Section 4.10. 


4.1.1 Sirius 

Like Bailey et al. (2015) we have used Sirius as a low polarisa¬ 
tion standard, and in fact combined the measurements reported 
there with additional observations. We now make its polarisation as 
5.5±1.7 ppm. From the discussion below in Section 4.3 it will be 
seen that this is entirely consistent with interstellar polarisation at 
southern declinations. At 2.6 pc distant and a V magnitude of -1.46 
Sirius is a particularly good unpolarised standard for instruments 
that can tolerate its brilliance. It should be pointed out though that 
the white dwarf companion, Sirius B, is currently separated from 
the primary by 10.2" (Burnham 1978a) and well outside the 6.7" 
HIPPI aperture. Care will need to be taken with larger instrument 
apertures or closer to periastron (which next occurs in 2044) when 
the separation is 3" (Burnham 1978a). 


4.2 Survey statistics 

In Table 6 the polarisation properties as a function of spectral type 
for both this work and the survey conducted by Bailey et al. (2010) 
with PlanetPol are listed for comparison. The PlanetPol instrument 
operated with a very broadband red filter covering the wavelength 
range from 590 nm to the detector cutoff at about 1000 nm (Hough 
et al. 2006) - redder than our measurements made with HIPPI. 
The empirical wavelength dependence of interstellar polarisation 
is given by Serkowski’s Law (Wilking et al. 1982): 

p(^/^...v) = exp((0. 1 - 1.86T„„) XnfAIA^a.)), (D 



p./p,. 

Figure 1. (a) How the polarisation measured, p, changes as the ratio of 
intrinsic to interstellar components, p,/p, . The range of values between the 
plotted maximum and minimum (dashed) are possible, but for a sufficiently 
large sample the mean (black) and median (grey) p will always greater than 
Pi for a non-zero p,. (b) As p,/p,- increases the probability of a random 
measurement of p for a star being greater than p, increases. 

where A is the wavelength examined and the wavelength 
of maximum polarisation. A typical value for is 550 nm 
(Serkowski et al. 1975), and for the effective wavelength of a GO 
star observed by HIPPI and PlanetPol this corresponds to a fac¬ 
tor 1.08 times greater polarisation for HIPPI. However, extremes in 
A,„ax can range from 340 to 1000 nm. The difficulty in measuring 
polarisation to the ppm level has so far prevented /l,„„j- from being 
determined within 100 pc. 

Before analysing Table 6 it is pertinent to briefly discuss the 
statistics of the degree of polarisation in such a survey as this. The 
degree of polarisation of a system, p, is the vector sum of intrinsic, 
p,, and interstellar, p,, components; from the Law of Cosines: 

P = ^P^'^i-'pj'^-^P;Prcos(7r^-'2^, (2) 

where a is the angle of interstellar polarisation relative to the direc¬ 
tion of the intrinsic component. As we do not know the orientation 
of either component of polarisation, a will have a random value 
between 0 and n. Figure 1 (a) shows how this effects polarisation 
measured as a function of the ratio of intrinsic to interstellar polari¬ 
sation. If the intrinsic component of polarisation is more than twice 
the interstellar component, then the polarisation measured will al¬ 
ways be greater than pj. However, even for relatively small values 
of p,, on the average we expect p to be greater than p/, even if a 
fraction of polarisations measured are less than the degree of inter¬ 
stellar polarisation - as shown in Figure 1 (b). 

From Table 6 the median polarisation is similar across the 
spectral classes A, F, G and K for both instruments. Tinbergen 
(1982) suspected the presence of variable intrinsic polarisation at 
the 10^“* level in stars with spectral type FO and later. This sup¬ 
position was not supported by the more sensitive observations of 
Bailey et al. (2010). Higher polarisations were indicated for spec¬ 
tral classes B and M, but they had very few stars of those types, 
and at larger average distances, and so could not ascribe that to any 
intrinsic polarisation of those stars. 

Adding our results to those obtained with PlanetPol makes it 
clear that B and M spectral classes do have higher polarisations; 
this is true even if we discount the Be stars. As is to be expected 
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from our selected V magnitude limit of 3.0 the B stars are on aver¬ 
age at larger distances than those of other spectral classes - more 
than twice the mean distance than for A stars - yet the median po¬ 
larisation is more than 4 times that of the A stars. We have fewer 
M stars, together with the PlanetPol only eight, hut three of them 
have a polarisation of ~500 ppm or greater. Additionally, the fur¬ 
thest K giant (K3III), 5 Sgr, also has a polarisation in excess of 500 
ppm. As will be discussed later, the interstellar contribution is not 
likely to be much more than 150 ppm even at 100 pc, and so a mea¬ 
surement of 500 ppm, in addition to the variability observed with 
respect to the observations of Tinbergen (1982), strongly argues for 
an intrinsic cause. 

Plotting the polarisation for all stars observed by us and Bailey 
et al. (2010) on a H-R Diagram (Figure 2) best illustrates the trends 
with spectral type. Even if one ignores the highest polarisations 
of the four Be stars, it can be seen that there is a stark contrast 
around spectral type AO, where earlier types have a greater degree 
of polarisation. The sharpness of the division between A and B stars 
suggests an intrinsic mechanism or mechanisms particular to this 
stellar type; this will be discussed in Section 4.9. 

On the giant branch (luminosity class III) we have very sparse 
data earlier than G5. For later types, and greater luminosities, there 
is a trend toward higher polarisations. In Section 4.10 we identify 
a number of K and M giants that we can be reasonably certain are 
intrinsically polarised and investigate the most likely mechanisms. 

4.3 Polarisation spatial distribution 

To investigate the spatial distribution of interstellar polarisation we 
have plotted PlanetPol and HIPPI survey stars’ polarisation degree 
and angle in equatorial co-ordinates (Figure 3) and projected on to 
Galactic co-ordinates (Figure 4). B type stars, M type stars and the 
K giant a Sgr have been neglected to reduce contamination from 
intrinsic polarisation. 

The PlanetPol survey (Bailey et al. 2010) found low interstel¬ 
lar polarisation toward the Galactic north pole, with higher polari¬ 
sations below 40° Galactic latitude toward the Galactic centre'’, and 
that in this range there was a tendency for alignment in the Galac¬ 
tic plane. Removing the B and M stars hasn’t changed that picture 
(Figure 4). The new results we present here align well with those of 
PlanetPol. The few data points we add at high Galactic latitudes are 
relatively low polarisation. Near the Galactic centre we add data at 
lower northern Galactic latitudes; the polarisation vectors for these 
data points between I = ~0-120° show an alignment consistent with 
the PlanetPol data, i.e. along the Galactic plane. Further south a dif¬ 
ferent picture emerges. 

Our survey has added significantly more data points in the 
Galactic south. Even if one allows for the wavelength difference 
between HIPPI using the g' filter and PlanetPol’s 600-800 nm win¬ 
dow, in general southern stars appear significantly more polarised 
than northern ones. That the degree of polarisation is greater in the 
Galactic south than north is consistent with the distribution of in¬ 
terstellar dust clouds seen around O, B, F and G stars by IRAS 
at 60 fjm (Dring et al. 1996). This might be explained by refer¬ 
ence to the Sun’s vertical displacement from the Galactic plane, Zq- 

® It should be pointed out that the spatial scale of measurements made by 
PlanetPol and HIPPI is many orders of magnitude less than the Galactic 
scale, and that in this case a maximum seen toward the Galactic centre is 
coincidental. The Solar System is currently located between the Orion Spur 
and the Sagittarius Main Arm. 


Many studies have shown that the Sun is slightly to the north of 
the Galactic plane compared to its assumed position in the Galactic 
co-ordinate system (Joshi 2007). Thus, by observing to the south, 
we are looking through more of the Galactic plane where we might 
expect interstellar grains to be more strongly aligned and impart 
greater polarisation. 

By plotting the polarisation degree per distance in Figure 5 
and with distance in Figure 6 we can refine the picture further. We 
denote the polarisation with distance as p/d, it is useful for com¬ 
paring the degree of polarisation of objects that might be separated 
by tens of pc. In Figures 5 and 6, the polarisation has been debiased 
as -sjp^ - cr^. Bailey et al. (2010) noted at right ascensions greater 
than n ha pi d relation of ~2xl0“^ pc“'. With the removal of the 
B and M stars from the map, there remain only 3 stars supporting 
this trend for declinations greater than 30°. One of these is Vega, 
which we can discount from considerations of interstellar polarisa¬ 
tion on the basis of its large debris disk. The remaining two stars 
are the K giants y Dra and k Lyr. From the discussion that follows 
in Section 4.10 it is likely that these stars are also intrinsically po¬ 
larised. Indeed, cross-reference between Figures 3 and 5 reveals 
the position angles of these three stars do not match those nearby. 
Given that they all lie near the Galactic equator, where there is oth¬ 
erwise a good measure of alignment, the case is strong for intrinsic 
polarisation. 

The large degree of polarisation observed for a Cen is prob¬ 
ably not real, and is discussed in detail in Section 3.1. From our 
survey, two further K giants, a Phe and a Tuc, can be identified 
as probably intrinsically polarised, as they have a p/d greater than 
any other stars within 100 pc. What remains is a region of low pid 
centred on 14 h right ascension, 1-35° declination - this region was 
mostly covered by the PlanetPol survey - indicating a relatively 
dust free volume. Within this region stars are polarised at ~2xl0“^ 
pc^F Further South and East p/d is greater, and somewhat patchy. 
The stars with the greatest p/d tend to lie within -15° and -30° 
declination or nearby; these stars are polarised at ~2xl0“* pc“*. 
Between the two regions are stars with intermediate polarisations. 
A linear fit to stars from our survey not suspected of being intrinsi¬ 
cally polarised gives 1.14x10“® pc“* with a coefficient of determi¬ 
nation, = 0.67. The highly polarised southern stars for the most 
part lie within 30 pc. Though we have fewer stars further away, 
these are mostly less polarised. This suggests that the higher inter¬ 
stellar polarisation of southern stars is imparted predominantly by 
nearby dust lying, from Figure 6, between 10 and 30 pc distant. Co¬ 
incidentally, the majority of estimates for Zq place the Sun 15 to 30 
pc above the Galactic plane (Joshi 2007), though the implication is 
most likely a local dust cloud at this distance. 

The polarisation of more distant stars in the interstellar 
medium has been shown to increase with distance as ~2xl0“® pc“' 
(Behr 1959). That nearby stars show less polarisation was demon¬ 
strated by Tinbergen (1982) and Leroy (1993b). Exactly how much 
less was determined by Bailey et al. (2010) for northern stars in the 
PlanetPol survey. Bailey et al. (2010) believed that the furthest stars 
in their survey showed high polarisations as a result of their prox¬ 
imity to the wall of the Local Bubble. As a result of the work pre¬ 
sented here, we now believe those stars to be polarised by intrinsic 
processes. The two furthest stars shown in the figures are Canopus 
and 6 Sco from our survey; they exhibit the greatest degree of po¬ 
larisation. This may be due to these stars being beyond the wall of 
the Local Bubble. Repeat observations of Canopus don’t show any 
variation (Bailey et al. 2015), which supports this. However, they 
are also the most luminous non-B stars in the survey, and to the 
best of our knowledge the polarisation of close bright giants at the 
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Spectral Type 


Figure 2. The results of both the HIPPI (green) and PlanetPol (red) surveys plotted on a H-R Diagram. The degree of polarisation detected is represented 
by the area of the bubble. The main sequence (V), sub-giant (IV), giant (III) and bright giant (II) branches are indicated. Two clear trends are observed: (1) 
increased polarisation of B stars, and (2) increased polarisation of late giant stars. 


ppm level has not previously been investigated (though it should be 
noted that early type supergiants do display intrinsic aperiodic vari¬ 
able polarisation arising as a result of asymmetric mass loss (Hayes 
1984, 1986)). The 3rd magnitude limit in V for this survey has re¬ 
sulted in fewer main sequence and sub-giant stars at distance than 
the PlanetPol survey with its limit of 5th magnitude. Without more 
data it isn’t possible to say definitively whether we’re probing the 
wall of the Local Bubble with our furthest stars. 


4.4 Debris disks 

Debris disks are circumstellar disks of dust around main-sequence 
stars. They are typically detected and characterised by excess in¬ 
frared emission above that of the stellar photospheric continuum 
emission. The peak wavelength of dust emission from debris disks 
depends primarily on the distance from the star (Matthews et al. 
2014). Broadly speaking, two typical temperature regimes are ob¬ 
served; warm dust (Tj 150 to 220 K) analogous to the Asteroid 
belt, and cold dust (T^ ~ 50 to 80 K) analogous to the Edgeworth- 
Kuiper belt (Morales et al. 2011). Edgeworth-Kuiper belt analogues 
are more easily detected by this method than closer disks owing to 
the greater contribution of the stellar photosphere to the total emis¬ 
sion at mid-infrared wavelengths. Recent far-infrared surveys with 
the PACS instrument (Poglitsch et al. 2010) on the Herschel Space 
Observatory (Pilbratt et al. 2010), where detection is already sen¬ 
sitivity limited, have found debris disks for up to 33% of A stars 
(Thureau et al. 2014), and 20% of FGK stars (Eiroa et al. 2013). 


The ratio of the luminosity of the disk to that of the star, Eir/E*, is 
indicative of the total starlight intercepted by the debris disk (Wy¬ 
att 2008), but by itself is not a reliable predictor of scattered light 
brightness (Schneider et al. 2014). The spectral energy distribution 
of a debris disk is also a function of the dust particle size distribu¬ 
tion and composition as well as its architecture. Dust particles in 
circumstellar disks polarise light by scattering and absorption pro¬ 
cesses, and sensitive polarimetry is potentially useful for removing 
degeneracies (Garcia & Gomez 2015; Schneider et al. 2014). Po¬ 
larisation seen by aperture polarimetry - where the aperture takes 
in the central star as well as the whole/a large portion of the disk - 
has been reported at levels of ~0.1 to 2% (Garcia & Gomez (2015) 
and references therein). 

As part of our survey we observed four main sequence ob¬ 
jects thought to be debris disk host stars, namely: Fomalhaut (HIP 
113368), I Cen (HIP 65109),/? TrA (HIP 77952) and y TrA (HIP 
74946). The PlanetPol survey also observed five debris disk sys¬ 
tems: Vega (BS 7001), Merak (BS 4295), y Oph (BS 6629), /? Leo 
(BS 4534) and a- CrB (BS 5793). Concurrently with the bright star 
survey we carried out another program where we observed debris 
disk systems. What can be learned from all of these observations 
has been substantially dealt with in another paper (Marshall et al. 
2015). In brief, comparison of the polarisation with measurements 
of the disks’ thermal emission demonstrate the capacity for polar¬ 
isation measurements to constrain the geometry and orientation of 
unresolved debris disks. A typical ratio of polarisation to thermal 
emission of between 5 to 50%, consistent with scattered light imag- 


MNRAS 000, H22 (2015) 
































10 


D. V. Cotton et al. 



Right Ascension, a(h) 

Figure 3. A to K stars within 100 pc from the HIPPI (green) and PlanetPol (red) surveys plotted on an equatorial co-ordinate system. Stellar positions are 
marked by dots, with a greyscale indicating distance from the Sun in pc. Vectors denote the degree and angle of the polarisation. 


100 ppm 



Figure 4. A to K stars within 100 pc from the HIPPI (green) and PlanetPol (red) surveys plotted on Galactic co-ordinates. Stellar positions are marked by dots, 
with a greyscale indicating distance from the Sun in pc. Vectors denote the degree and angle of the polaiisation. 
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Figure 5. A to K stars within 100 pc from the HIPPI (green) and PlanetPol (red) surveys plotted on an equatorial co-ordinate system. Stellar positions are 
marked by dots. Polarisation with distance, pjd, in ppm/pc is given by the colour scale to align with Figure 6. Brown rings denote debris disk stars. Stars 
labelled in red are K giants that may be intrinsically polarised. Stars labelled in black may not be representative of interstellar polarisation for other reasons 
(see the text). 



Figure 6. Polarisation with distance, pjd, for A to K stars within 100 pc from the HIPPI (green) and PlanetPol (red) surveys. Stellar positions are marked by 
dots. Polarisation with distance in ppm/pc is given by the colour scale to align with Figure 5. Brown rings denote debris disk stars. Stars labelled in red are 
K giants that may be intrinsically polarised. Stars labelled in black may not be representative of interstellar polarisation for other reasons (see the text). The 
dashed line represents a linear fit to the unannoted HIPPI stars. 
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Q 



Figure 7. The vector sum, p, of interstellar, pi, and intrinsic, p*, compo¬ 
nents of polarisation shown in the Q-U plane. The situation where the sum 
is less than the interstellai* polaiisation by itself is represented by the red 
arc; the situation where it is greater is represented by the blue arc. Given 
a random distribution more points will lie on the blue arc than the red and 
thus the average p will be greater than pj. The polarisation is equivalent to 
the median when p^ is perpendicular to pj, and then Equation 2 reduces to 


ing measurements. The majority of disks in the sample have po¬ 
larisation signals aligned approximately perpendicular to the disk 
major axis, indicative of scattering from the limh of small grains in 
those disks. No trend was found in the polarisation with either disk 
inclination, nor luminosity of the host star. The characterisation of 
the system j8 TrA as debris disk system based on a far infrared ex¬ 
cess derived from a single waveband IRAS measurement is found 
to be unreliable. 

What remains to be presented here is a comparison of polari¬ 
sation in these debris disk systems compared with those of the other 
main sequence stars in the survey, this is shown in Table 7. If the 
contribution from the disk is greater than the interstellar polarisa¬ 
tion then it is obvious that we would expect a higher polarisation on 
average from debris disk systems. Yet geometrical considerations 
related to the vector sum of intrinsic and interstellar components 
mean we should also expect a slightly higher polarisation from de¬ 
bris disk systems on average when the magnitude of contributions 
is similar or even smaller (as detailed in Figure 1). The situation is 
depicted graphically in Figure 7. 

The analysis of the PlanetPol survey data is presented here 
for the first time. In the PlanetPol data, debris disk systems show 
slightly higher median polarisations at shorter average distances. 
If one were to apply the pjd relation established by Bailey et al. 
(2010) for stars between 10 and 17 hours RA in the PlanetPol sur¬ 
vey, then one would expect 3.8 ppm at 19.0 pc. However as both the 
debris disk and main sequence samples contain stars with greater 
RA, a more cautious value of 5.1 ppm is achieved by scaling the 
pjd from the other class V stars. The relation described by Figure 
7 then gives us a median polarisation contribution from debris disks 
in the PlanetPol sample of 8.1 ppm. Given PlanetPol’s 5" aperture 


Table 7. Polarisation for debris disk and non-debris disk main sequence 
stars. 


Type 


HIPPI” 


PlanetPol* 


N 

Mean 

Median 

N 

Mean 

Median 



d(pc) 

P (ppm) 


d(pc) 

P (ppm) 

Debris Disk^ 

3 

27.4 

29.5 

5 

19.0 

9.6 

Other class V'^ 

6 

16.2 

18.1 

11 

33.0 

oq 

00 


a - This work, 
b - Bailey et al. (2010). 

c - If 7 TrA in HIPPI survey is excluded: 2, 12.9, 29.1. 
d - A-K stars only. Doesn’t include a Cen in HIPPI survey (see Section 
3.1.1); if included: 7, 14.1,27.5. 


wouldn’t often capture the full disk for such close objects this is in 
line with expectations. 

In the HIPPI data we have a higher polarisation for debris 
disks, but at a larger average distance. However y TrA has a tiny 
infrared excess and would not be expected to have much of a po¬ 
larisation signal, maybe 5 ppm at most (Marshall et al. 2015), and 
yet at 56 pc it is contributing most to the mean distance. The other 
two disk systems together have a median polarisation of 29.1 ppm 
at an mean distance of 12.9 pc, which given their median excess 
of Tir/T* = 18x10“* (Marshall et al. 2015) is more in line with 
expectations. However, the very small numbers of debris disk sys¬ 
tems (and also other main sequence stars) means that it is difficult to 
draw conclusions at this level, particularly in light of an uncertain 
contribution from the local interstellar medium (Section 4.3). 

Comparison of Figure 3 and Figure 5 reveals four of the eight 
debris disk systems to have polarisation angles closely aligned with 
their nearest neighbour stars in the survey. This naturally leads to 
two hypotheses: (1) that interstellar polarisation is swamping the 
contribution from the disk, or that (2) the local interstellar mag¬ 
netic field plays a role in the formation of the system. In this in¬ 
stance both of these may be discounted by consideration of the 
individual systems. The four systems with similar alignment are 
Fomalhaut, l Cen, Merak and y Oph. Taking the second hypothe¬ 
sis first, this would seem most probable if the systems were young 
and had formed near to where they are now. None of them is par¬ 
ticularly young though. Furthermore, the HIPPI aperture does not 
encompass the whole disk of Fomalhaut nor Merak. Meaning that 
the position angle we have measured does not correspond to the 
minor axis of the disk for either of these systems (Marshall et al. 
2015). 

The contribution of the interstellar polarisation to that mea¬ 
sured for the debris disk systems is harder to gauge without more 
precise data from nearby systems. Certainly, the data presented 
in Table 7 suggest that the interstellar contribution could be sig¬ 
nificant. A number of factors argue against interstellar polarisa¬ 
tion being dominant though. In the case of Fomalhaut the p/d is 
3.06 ppm pc“', higher than for any other survey star identified as 
having interstellar polarisation alone. For t Cen we have obtained 
data in different wavelength bands as part of another project (Mar¬ 
shall et al., unpub. data); preliminary analysis of these data gives 
a wavelength dependence inconsistent with interstellar polarisation 
(Serkowski 1973; Oudmaijer et al. 2001). For y Oph the polari¬ 
sation angle measured is very well aligned with the minor axis of 
the disk determined from imaging (Marshall et al. 2015), making 
the interstellar contribution very difficult to gauge. The PlanetPol 
measurements of Merak and p Leo are consistent with interstellar 
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polarisation -p Leo has a polarisation of only 2.3± 1.1 ppm and we 
have, on occasion, used it as a low polarisation standard - and it 
could be the characteristics of the disks in these two systems do not 
generate significant polarisation. For a fuller discussion of factors 
contributing to debris disk polarisation see Marshall et al. (2015). 

4.4.1 Altair 

Altair (BS 7557, a Aql), observed as part of the PlanetPol survey 
(Bailey et al. 2010), though not a classical debris disk system, has 
been identified as having an infrared excess through interferometric 
measurement of the star at near infrared wavelengths (Absil et al. 
2013). Based on the spectral slope of the excess it was predicted 
that the excess was attributable to scattered light rather than ther¬ 
mal emission. Given that small dust grains should be more effective 
scatterers at shorter wavelengths an appreciable polarisation would 
have been expected. The PlanetPol measurement is consistent with 
the contribution from the interstellar medium. From this we might 
infer that the grains responsible for scattering could well be too 
small to be effective polarisers, i.e. nanoscale dust grains as postu¬ 
lated by Su et al. (2013). 

4.4.2 p Pup 

Another non-traditional disk system is the bright giant p Pup (FlIP 
39757, F5II). It was identified as having a debris disk by Rhee et al. 
(2007), with a 60 ^m excess of 5±lxl0“^; it has an even more sig¬ 
nificant excess according to McDonald et al. (2012) of 160.7x 10“®. 
The degree of polarisation measured, 18.3 ppm, is consistent with 
what is expected from its infrared excess (Marshall et al. 2015). 
However, there are no imaging data available to gauge its geometry 
or confirm that the dust present is in the form of a traditional debris 
disk. The system is 18.3 pc distant, positioned on the border of the 
low and high polarisation regions in Figure 5, and so it is difficult to 
determine the relative contributions of the interstellar medium and 
any intrinsic component due to a disk. 

4.4.3 eSgr 

There is one other debris disk system in the survey not mentioned 
up to this point: e Sgr (HIP 90185). It has a high polarisation mea¬ 
surement of 162.9±4.4 ppm. It is an unusual system to be discussed 
in the context of debris disks on two counts: (1) the primary star 
is a 3.52 Mq B giant having spectral type B9.5III, and (2) it is a 
binary system with the 0.95 Mq companion (Hubrig et al. 2001) 
separated on a similar scale to the debris disk (Rodriguez & Zuck- 
erman 2012). The secondary orbits it at 106 AU, whilst the disk 
has been detected in IRAS 60 pm to have an excess of 4.5x 10“® by 
Rhee et al. (2007) and is presumed to be centred at 155 AU as a re¬ 
sult (Rodriguez & Zuckerman 2012). Even greater excess has been 
detected with Spitzer at 13 pm and 31 pm (Chen et al. 2014; Mit¬ 
tal et al. 2015), implying a closer disk. The 60 pm determination 
places the secondary well within the HIPPI aperture, but the disk 
on the edge of it. e Sgr was observed on 1/9/2014 when the seeing 
was particularly bad (~4"), so we probably have a large contribu¬ 
tion to the observed polarisation due to scattering from the debris 
disk. 

What makes this system particularly interesting is that in V 
band a simple calculation shows that the contribution of the sec¬ 
ondary to the light reaching the disk varies from ~0.1% of the total 
for the furthest part to ~3% for the closest. There is therefore an 


asymmetry in reflected light over the whole disk and thus in polar¬ 
isation as well. If the polarisation we see is a result of this asym¬ 
metry then the measured polarisation angle should be related to the 
position angle of the binary system. The position angle of the sys¬ 
tem was measured in March 1999 to be 142.3° (Hubrig et al. 2001). 
Our measured polarisation angle is 38.1±1.6°. When one considers 
that the secondary could have rotated by as much as ~11° (assum¬ 
ing a face-on system, a circular orbit and the formal limit on the 
uncertainty of the mass of component A), and that there must be 
contributions from interstellar polarisation, it is probable that we 
have measured a polarisation angle perpendicular to the position 
angle of the binary system. This is consistent with the stated hy¬ 
pothesis. If the system is inclined then the position angle of the 
binary system will not have swept through as great an angle, but 
we don’t know precisely the interstellar contribution nor the con¬ 
tribution resulting from any asymmetry associated with disk incli¬ 
nation, so the small difference between the position angle and the 
expected polarisation angle is not significant. More significant is 
that the measured infrared excess is only 11.0x10“® (Mittal et al. 
2015), and this shouldn’t be enough to generate this degree of po¬ 
larisation. However, the distance of the debris disk has been deter¬ 
mined from the IRAS measurement at 60 pm and the shorter wave¬ 
length measurements made with Spitzer imply either a closer disk 
which would produce a greater asymmetry, or a circumsecondary 
disk as has been suggested for the HD 142527 system (Rodigas 
et al. 2014). 

Based on the IRAS and Spitzer data we modelled the spec¬ 
tral energy distribution of the system using two blackbody com¬ 
ponents and obtained an fractional excess of 31x10“®, which is 
more consistent with what we would expect based on the polari¬ 
sation measured. The quality of this value is strongly dependent on 
the assumptions made regarding the temperature of the cold com¬ 
ponent and the stellar photosphere contribution. The stellar pho¬ 
tosphere was represented by a Castelli-Kurucz model (Castelli & 
Kurucz 2004) with an effective temperature of 10,000 K and a sur¬ 
face gravity, log g, of 4, and solar metallicity. This was scaled to the 
optical and near-infrared photometry from SIMBAD. In our model 
we fixed the temperature of the cold component to 85 K, such that 
its emission peaked at 60 pm - the wavelength of the longest re¬ 
ported flux density measurement. The warm component was fitted 
to the mid-infrared excess, which resulted in a temperature of 300 
K. The combined model overestimates the flux density at 31 pm 
compared to that reported by Chen et al. (2014). Due to the ab¬ 
sence of longer wavelength data constraining the peak of the cold 
emission, the total fractional excess is subject to large uncertainties. 

At this point it should be noted that e Sgr B was found as a 
result of a search of late-B stars showing high X-ray fluxes (Hubrig 
et al. 2001). Some X-ray binaries have been found to show variable 
polarisation (Clarke 2010) and so this is an alternative explanation 
for the polarisation observed. However, such detections have been 
rare, and generally for much stronger X-ray sources. Considering 
the magnitude of the flux asymmetry induced by the secondary, 
scattering from the debris disk seems the most likely mechanism. 
Indeed the X-ray activity might be an indication of dust accretion 
onto the secondary, which is another scenario suggested for HD 
142527 (Rodigas et al. 2014). 

Presuming the polarisation we see in this system is a result 
of scattering from the disk, then the behaviour with wavelength 
will be a function of the ratio from components A and B, as well 
as the properties of the grains in the disk. The polarisation will 
likely be greater at redder wavelengths than typical for a debris 
disk system. If the secondary is interior to the disk, as it orbits it 
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will race ahead of the disk, illuminating each section in succession 
as if a torch shinning upon it. Follow up observations of this object 
offer a unique opportunity to use polarisation to probe the radial 
homogeneity of the debris disk. Though it will take some time... 

4.5 Ap stars 

As an aside, from Figure 5 it is clear that e UMa (BS 4905) is 
more polarised than the stars around it. e UMa has spectral type 
AlIII-IVp and is the brightest chemically peculiar star in the sky. 
Ap stars typically exhibit high degrees of polarisation as a result of 
magnetic fields in the many hundreds of Gauss (Clarke 2010). Such 
fields generating broadband polarisation through, for example, dif¬ 
ferential saturation of Zeeman components as described by Leroy 
(1990). e UMa has a strong magnetic field, and a rotational axis at 
an angle to its magnetic axis that periodically brings its magnetic 
pole into line of sight; this is the likely cause of its high degree of 
polarisation. Flowever, recently it was proposed that the periodic¬ 
ity displayed (5.0887 d) by the star is a consequence of a 14.7 Mj 
companion orbiting at a distance of 0.055 AU (Sokolov 2008), and 
this could also be a contributor. 

There is one other Ap star in the PlanetPol survey: CVn. 

This star is known to have a regular time varying polarisation. Re¬ 
cent broadband measurements in a similar range to that of Planet¬ 
Pol have shown that the degree of linear polarisation varies between 
zero and more than 0.7% (Kochukhov & Wade 2010). It thus ap¬ 
pears a matter of chance that the recorded degree of polarisation for 
CVn by PlanetPol was just 8.8 ppm. 

4.6 Eclipsing binaries 

The work credited with bringing about the beginning of stellar po- 
larimetry in earnest is that of Chandrasekhar (1946)^. In this work 
he described how polarisation arising from free electron scattering 
at the limb of a star might be observed by using an eclipsing binary 
to break the symmetry of the stellar disc. It was shown that the limb 
would be most polarised with the azimuth of vibrations tangential 
to the limb, and that the polarisation would fall away quickly from 
there, becoming zero in the centre of the stellar disc. Usually the 
symmetry of the stellar disc renders this effect undetectable in aper¬ 
ture polarimetry, but an eclipsing binary breaks the symmetry. The 
magnitude of the limb polarisation varies depending on the type 
of star, with earlier types showing greater polarisation (Kostogryz 
& Berdyugina 2015). Measurements of polarisation across the disc 
of the Sun were first used to confirm the effect, and the Sun has 
a maximum limb polarisation of ~12% (Kostogryz & Berdyugina 
2015). 

There is one eclipsing binary in our survey, 6 Cap (HIP 
107556, A7III). This star also happens to be chemically pecu¬ 
liar. The G type secondary has an orbital period of 1.022789 days 
(Eggleton & Tokovinin 2008). It is noteworthy that the error asso¬ 
ciated with our measurement of 6 Cap is larger than any other in 
the survey. This could be the result of observing the system during 
a transit or transitioning into or out of secondary eclipse producing 
variable limb polarisation effects. However, the weather was vari¬ 
able on the night of the observation, and it is equally likely the large 
error is associated with reduced signal from patchy cloud. 

The star a CrB (BS 5793, AOV) observed in the PlanetPol 

^ The serendipitous discovery of interstellar polarisation was made in the 
process of searching for the Chandrasekhar Effect (Clarke 2010). 


survey is also an eclipsing binary of the Algol type (Eggleton 
& Tokovinin 2008), but doesn’t show any polarimetric behaviour 
from those observations of note. 


4.7 Close binaries 

Variable polarisation can occur in close binary systems as a re¬ 
sult of light scattered from material co-rotating in the system or 
from emission line effects associated with stellar winds or gaseous 
streams (Clarke 2010). A common property of many early-type 
close binary systems is the presence of a gaseous extrastellar enve¬ 
lope, possibly of proto-solar material (McLean 1980). The dynam¬ 
ics of a binary system render this envelope asymmetric, resulting 
in an intrinsic polarisation signal that varies with the binary phase 
according to the system geometry with respect to the observer and 
the polarigenic mechanism (McLean 1980). The requirement for 
extrastellar material resulted in a divide between evolved and un¬ 
evolved stars when investigated by Pfeiffer & Koch (1977). Having 
investigated systems with separations up to ~5 AU, they found bi¬ 
naries with unevolved stars, with few exceptions, do not show in¬ 
trinsic polarisation; whilst those with evolved stars are likely to be 
intrinsically polarised if the pair are separated by more than 10 so¬ 
lar radii (10 Rq). Pfeiffer & Koch (1977) hypothesised that the 10 
Rq divide resulted from insufficient material being present to gen¬ 
erate a large enough signal to be detectable by the instrumentation 
of the time (that being ~10“‘*). 

Within our survey there are a three B-type primaries with close 
companions, these being HIP 60718 (a Cru - as already mentioned 
in Section 3.1.2), HIP 100751 (a Pav) and HIP 65474 (Spica). All 
three show significant polarisations. 


4.7.1 Spica 

Spica is classified as BlIII-IV and has a B2V companion with an 
orbital period of just 4.0145898 days (Harrington et al. 2009). It 
is Beta Cephei-type variable star that varies in brightness with a 
0.1738-day period as a result of its outer layers pulsing. The bright¬ 
ness variability has been shown to be due to line-profile variability 
as a consequence of surface flows induced by fidal forces (Harring¬ 
ton et al. 2009). 

Spica was investigated by Pfeiffer & Koch (1977) and found 
not to vary polarimetrically - having an observed polarisation of 
0.03±0.01% it was designated as being unpolarised. However, as 
our initial measurements of Spica were intermediate of those of 
Tinbergen (1982) given in Table 2 - even if nominally within the 
error in ^ - we considered it an interesting target to follow up. Ta¬ 
ble 8 shows the three measurements we have made of Spica. From 
these it is clear that it is varying and thus intrinsically polarised. 


4.8 Be stars 

Be stars are defined as non-supergiants of B spectral type that have 
exhibited episodic Balmer line emission. The origin of the emission 
is attributed to the ejection of a gas circumstellar envelope (CSE) by 
the star (Domiciano de Souza et al. 2003). Be stars are rapid rotators 
that exhibit episodic mass and angular momentum losses as well as 
disc formation and dissipation (Domiciano de Souza et al. 2014a). 
The rotation rate is usually said to be 70-80% of the star’s critical 
velocity (Porter & Rivinius 2003). The stellar winds of Be stars are 
faster than those of ordinary B stars, particularly in earlier types. 
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Table 8 . Individual measurements of Spica. 


Date" 

UT* 

q (ppm) 

u (ppm) 

P (ppm) 

en 

24/5 

11:02 

-154.0 ± 2.5 

24.7 ± 2.5 

156.0 ± 2.5 

85.4 ± 0.9 

29/6 

10:32 

-185.1 ± 2.9 

30.6 ± 2.8 

187.6 ± 2.9 

85.3 ± 0.9 

29/6 

12:49 

-188.1 ± 3.6 

57.2 ± 3.9 

196.6 ± 3.8 

81.6 ± 1.1 


a - All dates are 2015. 

b - The time given as hh:mm and is that corresponding to the beginning of 
the measurement from the third telescope position angle in the sequence. 
The target was reacquired for each position angle leading to some minor 
variation in the timing. There was also a short pause between the middle 
two measurements of the third observation owing to passing cloud. Each 
measurement represents 24 mins of data plus acquisition time. 


The winds are asymmetric, characterised by fast tenuous winds at 
the poles and stronger slower winds at equatorial latitudes. 

Such stars exhibit varying polarisation on both short and long 
timescales. On minute to hour timescales fluctuating polarisation is 
attributed to ejection events, where a ‘blob’ of material is formed 
at the stellar surface, distorting its shape, which co-rotates with the 
star with the Keplerian velocity of the inner decretion disc (Carciofi 
et al. 2007). The extra equatorial material adds to the already dis¬ 
torted shape of the star caused by its sub-critical rotation (which 
is also responsible for gravity darkening at equatorial latitudes) 
(Domiciano de Souza et al. 2003). On longer time scales a circum- 
stellar decretion disc of ionized gas is built up around the star, the 
presence of which also produces an infrared excess and a non-zero 
(usually dominant) polarisation signal according to its density and 
geometry - the polarisation angle is aligned with the rotational axis 
of the star (Carciofi et al. 2007). During a quiescent phase the disc 
is dissipated by radiation pressure and partial reaccretion onto the 
star (Carciofi et al. 2012). Be stars are known to be variable on time 
scales ranging up to decades. A quiescent Be star might appear as 
an ordinary B star before becoming active within just a couple of 
days (Peters 1986; Barnsley & Steele 2013). A subclass of these 
objects is the ^sheu star, where the emission lines are depressed as 
a result of viewing geometry, i.e. they star is equator-on to the ob¬ 
server (Saad et al. 2012). However, 'Bsheii behaviour is often seen in¬ 
termediate of that of B and Be behaviour, suggesting the CSE may 
sometimes develop, at least initially, with an imperfect alignment. 
In such systems polarimetry is a relatively more sensitive probe of 
disc development. 

The definitive statistical study of Be stars and their polarisa¬ 
tion is that of Yudin (2001); Table 9 is a reproduction of the data he 
collated for the Be stars in our survey. 

Our measurements of a Col (HIP 26634), a Ara (HIP 85792) 
and 7j Cen (HIP 71352) in the g' agree in position angle with those 
collated by Yudin (2001). For a Ara Meilland et al. (2007) re¬ 
solved the equatorial disc using VLTI/AMBER data and calculated 
models with orientation in agreement with the position angle given 
by Yudin (2001), after subtraction of an interstellar component of 
0.15% at 30'’ to give 166°. Given the discussion in Section 4.3, this 
determination of the degree of interstellar polarisation seems high. 
Meilland et al. (2012) conducted a spectro-interferometric survey 
of Be stars which included both a Ara and a Col. The position an¬ 
gle for the rotational axis of these stars was obtained by means of 
a axisymmetric kinematic model. The determined position angles 
were 88±2° and 10° for a Ara and a Col respectively. This places 
our determined position angle for the circumstellar disc of a Ara 
between that of Meilland et al.’s 2007 and 2012 determinations. 


Table 9. Observed polarisation for Be stars from Yudin (2001) compared to 
our measurements. 


Stai' 

Yudin (2001)" 

This Work 

HIP 

P(%) 

en 

P(%) 

en 

7588 

0.04 

136 

0.21518 ±0.00038* 

31.2 ±0.1 

26634 

0.15 ± 0.05 

109 

0.07220 ± 0.00064 

101.1 ±0.5 

71352 

0.61 ± 0.04 

174 

0.53206 ± 0.00038 

172.0 ± 0.0 

85792 

0.37 ± 0.07 

174 

0.62400 ± 0.00034 

171.8 ±0.0 


a - Note that the tabulated values are the observed values, and that for 
some objects Yudin (2001) also calculated the intrinsic and interstellar 
components. 

b - Tabulated value for this work is the average of two measurements. 

and our measurement for a Col in closer agreement with Meilland 
et al. (2012) than that of Yudin (2001) is. 

Our measurement for the degree of polarisation of a Col is half 
the tabulated value indicating a smaller circumstellar disc, whereas 
for rj Cen we record almost double the polarisation level, indicating 
a thicker circumstellar disc. For a Ara we record a similar level 
of polarisation. There is no agreement for Achemar (HIP 7588) 
though, and as we observed it twice it bears closer scrutiny. 

4.8.1 Achemar 

As one of the closest brightest Be stars Achemar has been ex¬ 
tensively studied from the 1970s onward. During that time it has 
been in and out of emission - gaining and losing its disc. Dur¬ 
ing the most recent period quiescent period Domiciano de Souza 
et al. (2014b) took advantage of a negligible disc in using inter¬ 
ferometry to measure the position angle of Achernar’s rotational 
axis as 216.9±0.4°. They also measured the rotational flattening of 
Achemar as R^^IRp = 1.352±0.260 making it one of the flattest 
fast-rotating stars. When the star is in an emission or shell state it 
will have a significant disc that would be expected to lie with its ma¬ 
jor axis perpendicular to the star’s rotational axis. 36.9±0.4° is thus 
the angle one would expect from a polarimetric measurement of the 
star in a disc bearing phase before any other considerations. His¬ 
toric measurements of the polarisation of Achemar together with 
our new measurements are given in Table 10. 

McDavid (2005) interpreted the results of Schroder (1976) as 
being consistent with measuring the rotationally oblate photosphere 
of Achemar during a phase where the circumstellar disc was com¬ 
pletely absent. However, the degree of polarisation recorded, even 
allowing for uncertainty, is more than double what one would ex¬ 
pect from the calculations of Sonneborn (1982) given Achernar’s 
established inclination angle of ~60° (Domiciano de Souza et al. 
2014b) to ~65° (Carciofi et al. 2007). Alternative explanations are 
that this is the interstellar component of the observed polarisation, 
or that in the absence of the CSE it might be due to Achemar’s po¬ 
lar wind (Stee & Meilland 2009) or a combination of one or more 
of these. 

Excepting that of Schroder (1976) the data in Table 10 fall 
into two distinct categories: those where there was a significant disc 
present, and those where there was, perhaps, only a tenuous disc. 
The distinction being approximately an order of magnitude in p. 
The earlier position angle determinations are consistent with that 
expected from the interferometry measurements made by Domi¬ 
ciano de Souza et al. (2014b). However, our measurements, and the 
averages of those made by Carciofi et al. (2007) are more precise 
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Table 10. The polarisation of Achemar over time. 


Date 


P(%) 



en 


Band 

Reference 

31/1/1968 

0.03 

± 

0.01 

26 

± 

10 

V 

Serkowski (1970) 

12/12/1968 

0.02 

± 

0.01 

46 

± 

14 

V 

Serkowski (1970) 

1969-70" 

0.011 

± 

0.004 

136.2 

± 

10.4 

V 

Schroder (1976) 

1977-78" 

0.001 

± 

0.005 




400-700 nm 

Tinbergen (1979) 

1/9/1995'’ 

0.14 

± 

0.04 

41 

± 

8 

V 

McDavid (2005) 

1/9/1995'’ 

0.11 

± 

0.05 

45 

± 

13 

B 

McDavid (2005) 

7/2006" 

0.159 

± 

0.003 

31.7 

± 

0.6 

V 

Carciofi et al. (2007) 

9-11/2006" 

0.130 

± 

0.001 

32.6 

± 

0.3 

V 

Carciofi et al. (2007) 

9-11/2009""' 

0.011 

± 

0.013 




B 

Domiciano de Souza et al. (2014b) 

29/6/2011"' 

0.016 

± 

0.038 




B 

Domiciano de Souza et al. (2014b) 

9-10/2011""' 

0.020 

± 

0.015 




B 

Domiciano de Souza et al. (2014b) 

9-11/2011""' 

0.015 

± 

0.009 




V 

Domiciano de Souza et al. (2014b) 

1/7/2012"' 

0.017 

± 

0.015 




V 

Domiciano de Souza et al. (2014b) 

21 /11/2012"' 

0.035 

± 

0.050 




V 

Domiciano de Souza et al. (2014b) 

2/9/2014 

0.26984 

± 

0.00045 

31.5 

± 

0.1 

g' 

This work. 

24/5/2015 

0.16052 

± 

0.00062 

31.8 

± 

0.2 

g' 

This work. 


a - Exact dates of observations not reported. 

b - Measurements for U, R and I bands also reported in McDavid (2005) are the same within error as the B and V band measurements. 

c - For the sake of brevity we have averaged similar measurements reported by Carciofi et al. (2007) and Domiciano de Souza et al. (2014b). In the case of the 
former the measurements were weighted by the number of wave-plate positions indicated. Where p/cr > 4 debiasing has been caiTied out as 
otherwise the median is given for p, as recommended by Clarke (2010). 
d - Position angles not reported. 


and disagree by ~5°. Carciofi et al. (2007) observed the polarisa¬ 
tion to change in both degree and angle on a timescales of days 
to weeks, as well as on hour long timescales. They interpreted the 
change in position angle as likely associated with a departure from 
a simple 2D geometry for the disc. Changes in position angle of the 
order of 5°, and in the degree of polarisation of ~0.02% were seen 
to occur on timescales of less than an hour. 

Here our values represent observations made over a little less 
than half an hour each, so it’s reasonable to assume we’ve seen a 
distorted disc geometry due to the addition of transitory hlobs. Our 
measurement of 2nd September 2014 represents the highest polari¬ 
sation ever measured for Achemar (the previous highest being that 
of Carciofi et al. (2007) on 14th July 2006 of 0.188±0.009%) and 
thus the greatest recorded build-up of circumstellar material. By the 
24th of May 2015 the polarisation level had dropped by a third, but 
was still, historically speaking, very high. 

From Table 10 our measurements are far more precise than 
those recorded previously. HIPPI therefore has a previously unseen 
capacity to probe the growth and dynamics of circumstellar discs 
around Be stars. The growth and dynamics of circumstellar discs 
are still poorly understood (Carciofi et al. 2007), high precision po- 
larimeters like HIPPI clearly offer an enhanced capacity to study 
this phenomenon. 


4.9 B stars 

Setting aside the Be stars Table 6 makes it clear that B stars are 
more polarised than A, F, G, and K type stars. The starkness of 
the transition from A stars is shown up particularly well in Figure 
2; it also appears that early B stars are more polarised than later 
B stars. After excluding the Be stars, the close binaries and the 
unusual debris disk system e Sgr, there are five stars of B spectral 
type remaining. Although the case for y Crv and a Mus (which 
is 96.7 pc distant) is marginal, they all have polarisations that are 


higher than would be expected from only interstellar polarisation 
when considering their positions in Figures 5 and 6. 

Early type stars can be expected to have higher magnetic fields 
in general as a consequence of their younger age. Some Bp (Briquet 
2014) and Be stars (Hubrig et al. 2013) have been shown to have 
significant magnetic fields (though the fields in Be stars are weak 
- less than '-100 Gauss). However, no magnetic fields have been 
detected in ‘normal’ B stars within a margin of error of 13 Gauss 
(Shorlin et al. 2002). Consequently we don’t expect these stars to 
show polarisation from differential saturation of Zeeman compo¬ 
nents or similar mechanisms. 

One possibility to account for the high degree of polarisation is 
that these might be debris disk host stars. Infrared excess associated 
with debris disks is preferentially found in young systems. How¬ 
ever, Padgett et al. (2012) have recently, from WISE data, given the 
portion of B stars within 120 pc hosting debris disks as 12.1±2.0%. 
Similarly Schultz & Wiemer (1975) found an infrared excess - in¬ 
dicative of a debris disk system - in 20 of 350 O and B stars. Given 
that we have already identified one such system in our survey, the 
probability of many more is low. 

A feature of B type stars is that many of them have a high ro¬ 
tational velocity. From Table 2.1 it can be seen that a Gru, jj Lib 
and cr Sgr all have quite high rotational velocities. The spread in 
distances of these objects is quite large, and so we’ve plotted pjd 
against rotational velocity in Figure 8 to attempt to reduce the ef¬ 
fects of interstellar polarisation. Figure 8 shows B stars with higher 
rotational velocities tend to have higher polarisations. It appears 
that there is an exponential increase in pjd with V sin i. There are 
no exceptions to this trend in our data, but there are two stars sur¬ 
veyed by PlanetPol that don’t fit. 

T Her has a high pjd, especially amongst stars from the Plan¬ 
etPol survey. However, it lies at 94.6 pc distance, where we have 
very few stars, and the polarisation of the interstellar medium is not 
well constrained. The degree of polarisation we measured could be 
indicative of the edge of the Local Bubble - as originally suggested 
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Figure 8. Polarisation with distance, pjd as a function of rotational velocity 
for B stars without other identified polaiigenic mechanisms. For the stars in 
our survey there is a clear trend in increased polarisation with rotational 
velocity, whilst for the stars in the PlanetPol survey two objects don’t fit 
the trend: r Her and 6 Cyg. Parameters for the PlanetPol stars come from 
Bailey et al. (2010). 


for this object by Bailey et al. (2010). So, the polarisation measured 
could be interstellar in origin, but the lack of other objects at similar 
distances in the same part of the sky makes it hard to draw conclu¬ 
sions. The polarisation of 6 Cyg is not so easily explained away 
as being interstellar. It’s a multiple system but the B component 
(FIV) is sufficiently separated that we would not expect this to be 
the polarigenic mechanism in this case. It could be that S Cyg has 
a high rotational velocity, but a low inclination giving a low value 
for V sin i, but this is purely speculation. 

There are two candidates for a polarigenic mechanism related 
to rotational velocity. We begin by investigating that proposed by 
Bailey et al. (2010) for Regulus (BS 3982) - the star with the high¬ 
est V sin i in Figure 8. 


4.9.1 Rotational oblateness 

Based on the work of Chandrasekhar (1946) on revealing limb po¬ 
larisation with an eclipsing binary used to break the symmetry of 
the disc, Ohman (1946) reasoned that a rotationally flattened star 
must show polarisation at all times. Detailed modelling of the ex¬ 
pected behaviour has since quantified the effect expected, with the 
most recent predictions relevant to broad band polarimetry being 
made by Sonneborn (1982). The magnitude of the polarisation ex¬ 
pected is now believed to be much lower than originally suggested, 
and in almost 70 years of searching, polarisation induced by elec¬ 
tron scattering in a rotationally flattened atmosphere has not yet 
been confirmed. 

A degree of polarisation of 500 ppm parallel to the rotational 
axis is expected at 450 nm for a BO star inclined at 90° and rotating 
at 95% of its critical velocity. This amount is greatly reduced for 
inclination angles more face on. It also decreases with increasing 
spectral type, with the maximum in the visual also shifting redder. 
The progression is dramatic with a B5 star expected to show at most 
50 ppm at 700 nm (Sonneborn 1982). 

As already mentioned, Bailey et al. (2010) have proposed this 
mechanism for the B6 star Regulus, which rotates at 86% of its crit¬ 
ical velocity (McAlister et al. 2005). Regulus is close by in a part 


of the sky with low interstellar polarisation and registered 36.8± 1.6 
ppm at an effective wavelength of 741.2 nm with PlanetPol. We be¬ 
lieve Regulus still to be the best candidate for demonstrating polar¬ 
isation due to rotational oblateness, but unless the effect has been 
significantly underestimated it cannot explain the high polarisation 
of the other B-type stars seen here. 

We see larger polarisations in g' than were seen by PlanetPol 
in its redder pass-band (~575-1025 nm). There is no trend with 
spectral type in Figure 8, and our stars are rotating much slower 
than Regulus. Even accounting for a more significant interstellar 
contribution to the observed polarisation, the excess observed is far 
too large to be attributed to rotational oblateness. 


4.9.2 B stars as Be stars 

Another possibility is that the B stars that show significant polari¬ 
sation are actually Be stars outside of an emission phase, or that the 
polarisation we see is the result of a kind of sub-Be behaviour, or 
perhaps evidence for a polar wind in the absence of a CSE. In par¬ 
ticular if a small amount of ejected gas is present at the equatorial 
regions this could produce the observed polarisation. This hypoth¬ 
esis is advanced on account of the rotational velocity trend (Figure 
8) as well as the sharp division between A and B type stars in the 
H-R diagram (Figure 2). The statistical analysis of Yudin (2001) 
shows a triangular distribution for Be star polarisation with V sin i 
that peaks at 200-250 km/s, we have far fewer B stars, but the dis¬ 
tribution is not inconsistent with that picture. It has been suggested 
that faster rotators are more likely to have polar winds and as a 
result a circumstellar disc (Stee & Meilland 2009). 

A general description of Be behaviour and the transition from 
B to Bshell to Be has already been given in Section 4.8. In early type 
Be stars the strong winds are enough to drive material away from 
the star creating the CSE. That the winds of B stars are weaker 
might be a problem for this hypothesis, but later type Be stars have 
weaker winds and still show Be behaviour. Other mechanisms like 
non-radial pulsations, magnetic activity and binarity have been ad¬ 
vanced to explain CSE development in later Be stars (Saad et al. 
2012; Stee & Meilland 2009). 

It is worth keeping in mind that the levels of polarisation we 
are seeking to explain are only ~100 ppm - the same level recorded 
by Schroder (1976) for Achemar attributed to a disc-free state. This 
is more than three orders of magnitude lower than the level mod¬ 
elled by Wood et al. (1997) for the Be star rj Tau for instance. 
Naively scaling the amount of infrared excess they modelled for 
their thin disc model in correspondence with the degree of polari¬ 
sation we see here would result in something very difficult to detect 
given that infrared excess is attributable to gas and dust. Indeed 
from Table 2.1 it can be seen that the excess from B stars is not too 
different from that of Be stars in the survey. So, it may well be that 
polarisation is now a more sensitive probe of close circumstellar 
gas than infrared excess. It is certainly more sensitive than emission 
line measurements, as when there is little gas there is no emission. 
Furthermore polarisation is more sensitive to the inner part of any 
disc, whereas infrared is a better probe of the outer parts (Wood 
et al. 1997). We could be seeing a tenuous proto-CSE, and while 
a weak wind might not prevent a gas disc forming, it would likely 
prevent it developing as readily. 

Recently the issue of transition from B to Be star was looked 
into by Barnsley & Steele (2013). They examined Ha variability 
in a representative sample of Be stars. They formed a number of 
conclusions: (1) that the full phase of transition between B and Be 
star probably occurs over centuries, (2) that stars with earlier spec- 
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Table 11. Repeat observations of cr Sgr. 


Date 

P (ppm) 

PA(“) 

1/9/2014 

175.1 ± 3.8 

130.5 ± 1.3 

22/5/2015 

164.5 ± 5.9 

129.0 ±2.1 


tral type and higher V sin i show greater variability in Ho- emission, 
and (3) that for stars with smaller V sin i or later spectral type vari¬ 
ability was more likely on longer (years) timescales. What we take 
away from this is that stars with a higher rotational velocity or an 
earlier spectral type are more likely to show the emission character¬ 
istics of classical Be stars. But that any B star might become a Be 
star at any point in time, and that for slower later stars the behaviour 
is more subtle. 

It stands to reason then that B stars with higher rotational 
velocities are also more likely to form the gas shells needed for 
emission behaviour, even if insufficient gas is actually stockpiled 
to generate the emission. If this is the case we’d expect these shells 
to form at equatorial latitudes and result in polarisation perpendic¬ 
ular to the rotational axis just like in Be stars Porter & Rivinius 
(2003); Saad et al. (2012); Domiciano de Souza et al. (2014b). Fur¬ 
thermore, whilst variations in degree of polarisation are expected, 
variations in the polarisation angle should be small - as they are 
with Achemar for instance. This hypothesis has the advantage of 
explaining the trend we see with V sin i but also allowing for the 
exceptions of t Her and S Cyg. 

A pair of observations of cr Sgr (HIP 92855, B2V) shows 
the expected Be-like behaviour, albeit at a marginal level (Table 
11). Our observations therefore represent evidence for sub-Be be¬ 
haviour in ‘ordinary’ B stars. Incidentally, the inclination of cr Sgr 
has been estimated as 42±4°, putting its rotational velocity at 60- 
78% of critical (Hutchings et al. 1979)* - in line with typical Be 
stars. Higher than normal X-ray emission for cr Sgr has also been 
seen hy XMM Newton (Oskinova 2012), which according to Hubrig 
et al. (2013) is a behaviour associated with the Be phenomenon. 

After cr Sgr the next most polarised ‘ordinary’ B star in the 
survey is p Lib (HIP 74785, B8Vn). Notably it has the highest in¬ 
frared excess amongst these stars (Table 2.1). The inclination of 
P Lib has been estimated as 59j;g'° (i.e. 45-64% of critical veloc¬ 
ity) by Hutchings et al. (1979). The case for p Lib being a Be 
star is strengthened by its listing in SIMBAD as a variable star - 
rare/occasional low-level photometric variability being associated 
with fading in late Be stars (Hubert & Floquet 1998) - if it is vari¬ 
able it must be so on long timescales as it has been recommended as 
a photometric standard by a variety of later sources (e.g. Adelman 
(2001), references within Baade (1989)). Additionally McLaughlin 
(1932a) made a rather extraordinary aside in examining plates of p 
Lib from 1912 to 1915 (to test a claim that it was a spectroscopic 
binary^), saying, “Only the Balmer series of hydrogen is at once 
evident, and these lines are so broad and dilTuse that measurements 
are difficult, even on over exposed plates. On some spectrograms 
H^ and H^,- appear triple, as if weak emission were present.”'^ 


° Note the star name is incorrectly given as 6 Sgr in that work. 

® Even today it is uncertain whether/? Lib is a binary - Roberts et al. (2007) 
have identified a second star ~2.1" away, but they indicate follow-up obser¬ 
vations are needed to confirm it isn’t a background star. If it is a companion 
its spectral type is estimated as M2V. 

This statement gains greater significance when one realises that 
McLaughlin authored two other papers on Be star spectroscopy in the same 


4.10 Late giants 

The survey also contains a significant fraction of late-type giants. 
Two of these have significant previously measured polarisations: 
HIP 89931 (d Sgr, K3III) at the 10-“ level (Tinbergen 1982) and the 
SRB type pulsating variable HIP 112122 (fi Gru, M5III) also at the 
10-“ level (Heiles 2000). Both of these stars also show polarisations 
at that level here. As discussed earlier in Section 4.1 our measure¬ 
ments are significantly different to the earlier ones, indicating that 
the polarisation is variable and therefore intrinsic to the star. In ad¬ 
dition HIP 79593 (d Oph) has a similarly high polarisation. At 61.3 
pc distance a polarisation of 131.1 ppm for HIP 110130 (a Tuc, 
K3III) might suggest that it is also intrinsically polarised. When 
one considers that in the PlanetPol survey that the MO giants a Vul 
(BS 7405) and y Sge (BS 7635) recorded 1321.4 ppm and 199.5 
ppm respectively, this amounts to four out of eight luminosity class 
III M giants from the two surveys that are probably intrinsically 
polarised, and very probably two K type class III giants that are as 
well, based purely on the degree of polarisation. In addition, careful 
examination of the spatial distribution and p/d relations in Section 
4.3 suggests that y Dra (BS 6705, K5III) and k Lyr (BS 6872, K2III) 
may also belong to this group. Likewise a Phe (BS 2081, K0.5IIIb) 
has a polarisation/distance of 3.17 ppm/pc, higher than any other 
star we consider to show interstellar polarisation only. 

From Figure 2 we see that the degree of polarisation tends to 
increase with spectral type and also with luminosity in these stars. 
Dyck & Jennings (1971) examined polarisation in 55 K and M gi¬ 
ants and supergiants. They found that almost all supergiants dis¬ 
played significant polarisation and that class III giants displayed an 
increasing likelihood of being polarised with spectral type from M3 
later. However, the sensitivity of their study was no better than IQ-* 
to 10-“. Given the results presented here, it now appears that many 
warmer giants are also intrinsically polarised, just to lower levels. 
Indeed, the study by Dyck & Jennings (1971) includes d Oph, a 
Vul and y Sge determining no intrinsic polarisation in each case". 

Dyck & Jennings (1971) studied the wavelength dependence 
of the polarisation in these stars, finding that it increased as /f-' 
before tailing off at blue or violet wavelengths, where the exact 
peak varied from star to star. For an MO spectral type the effective 
wavelength for HIPPI is 2.094 prcT^ whereas for PlanetPol it was 
1.264 ytim-'. This means that HIPPI is 1.66 times as sensitive to this 
phenomenon as PlanetPol was. Indeed, if one applies this factor to 
the median polarisation given for M stars by PlanetPol in Table 6, 
the value is closer to that obtained by our survey. 

4.10.1 Circumstellar dust 

Dyck & Jennings (1971) also noted that all the stars in their survey 
with significant polarisations were variable stars and that they all 
had an infrared excess at 11 pm (which is associated with TiO dust, 
TiO being one of the first molecules formed in a cooling stellar en¬ 
vironment). Some variability was common, and was, in fact, one 
of the criteria used to assign polarisation as intrinsic. The combina¬ 
tion of these properties resulted in the cause of the polarisation to be 
attributed as scattering or absorption from a circumstellar shell of 
solid condensates (Dyck & Jennings 1971; Jennings & Dyck 1972). 

year (McLaughlin 1932b; Curtiss & McLaughlin 1932), and was thus well 
practised in identifying emission characteristics in the photographic plates 
of the time. 

" Dyck & Jennings (1971) also measured another three stars from the 
PlanetPol survey without significant polarisations. 


MNRAS 000, ffi22 (2015) 





Polarisation of southern bright stars 19 


1000 . 



Fractional Luminosity, F^p^/Fi, (x10‘°) 


500 - 


4 12 / 5/2014 
♦ 24 / 5/2015 


400 - 


Q . 300 - 


200 - 


100 - 


5 


0.5 1.0 1.5 

Inverse Wavelength, A' 


(ijm-'') 


Figure 9. Polarisation with infrared excess in giants later than G5. Note 
that HIPPI is 1.66 times as sensitive as PlanetPol at detecting polarisation 
caused by dust in red giants, but that it also sampled a region with a higher 
interstellar polarisation. Fractional luminosities come from McDonald et al. 
(2012), other parameters for the PlanetPol stars come from Bailey et al. 
( 2010 ). 


In this picture, dust forms in the circumstellar environment around 
giant stars through a two-step process: pulsations in the outer enve¬ 
lope elevate molecular gas to the surface, whence the stellar wind 
blows it radially outward to condense as grains in a cooler environ¬ 
ment (Morris 1987). Dust is associated with IR excess and many 
red giants have significant IR excess (c.f McDonald et al. (2012)). 
For polarisation to be induced from the dust requires that the cir¬ 
cumstellar envelope be asymmetric with respect to the observer. 
There are a number of potential causes of the asymmetry which 
are debated. It has been shown that the wavelength dependence of 
polarisation due to circumstellar dust depends on particle size and 
composition (McCall & Hough 1980; Raveendran 1991) as well 
as the geometry of the system. Silicate or carbonaceous grains are 
most common, however polarisation may also result from scatter¬ 
ing from inner shells of condensates able to form in hotter regions 
such as corundum and iron-poor silicates which are transparent at 
NIR wavelengths (Ireland et al. 2005). 

In all, there are 24 giants, a bright giant and a supergiant later 
than G5 in both surveys. Of these 18 appear in the infrared excess 
catalogue of McDonald et al. (2012). All but two have an excess 
greater than 1x10“* in fractional luminosity, with a range encom¬ 
passing 11 /jm, peaking around 12 yum or shorter. We plot polarisa¬ 
tion with infrared excess in Figure 9. This reveals that the degree of 
polarisation is not well correlated with infrared excess. This does 
not necessarily mean that dust isn’t a prominent scatterer in the 
outer regions of these stars, just that the amount of dust does not 
account for the degree of polarisation seen. Geometrical considera¬ 
tions may be dominant, which might indicate that the stellar winds 
of later giants are less symmetrical, and consequently produce in¬ 
creasingly asymmetric dust-discs. In the model of Johnson & Jones 
(1991) an asymmetric wind develops with a star’s entry onto the 
asymptotic giant branch (AGB). Our observations here might in¬ 
dicate that it occurs sooner. Alternatively, a cooler atmosphere al¬ 
lows dust to form closer to the photosphere, which would produce 
a larger polarisation signal due to an increased scattering density. 


Figure 10. polarisation measurements for <5 Oph with inverse wavelength. 
The dashed guide-line is drawn from the origin through the first data point. 
The second data point lies close to the guide-line, showing a relation¬ 
ship. 


Table 12. Individual measurements of 6 Oph. 


Filter" 

Date 

q (ppm) 

u (ppm) 

P (ppm) 

0 (°) 

g' 

12/5/14 

-15.4 ± 

6.0 -487.2 ± 

4.9 487.5 ± 5.5 

134.1 ± 0.7 

F 

24/5/15 

-44.5 ± 

9.1 -391.2 ± 

9.2 

393.7 ± 9.1 

131.8 ± 1.3 

g' 

24/5/15 

-84.8 ± 

7.2 -473.0 ± 

7.1 

480.6 ± 7.2 

129.9 ±0.9 

425SP 24/5/15 

-109.4 ± 28.3 -433.1 ±28.2 

466.7 ± 28.2 

127.9 ± 3.6 


a - The effective wavelength for each filter in //m ' is: 425SP: 2.502, g ': 
2.094, /: 1.654. 

4.10.2 S Oph and other mechanisms 

To investigate further we followed up our initial May 2014 obser¬ 
vation of the suspected variable (Percy & Shepherd 1992) d Oph 
with observations made in three different wavelength bands in May 
2015. This object is particularly interesting because it features in 
the work of both Tsuji (2008) on MOLspheres (extended molecu¬ 
lar spheres thought to exist beyond the photosphere) and the work 
of Ryde et al. (2015) which largely refutes the need for a MOL- 
sphere to explain spechal features such as metallic emission lines, 
OH lines and HF lines at 12 pm. Ryde et al. (2015) claimed that all 
of their sample, including 6 Oph, were dust free, and used this as 
an argument against certain absorption features being produced in 
the MOLsphere against a continuum formed from an alumina dust 
shell. 

The measurements of 6 Oph taken a year apart produced the 
same degree of polarisation within the error of the measurement as 
shown in Table 12. However, the polarisation in q was substantially 
different: -84.8±7.2 ppm for the later measurement compared to 
-15.4±6.0 ppm for the earlier one. The polarisation of 8 Oph is 
therefore varying, even if the large, nearly static, value for u masks 
the variation if just looking at the polarisation degree or angle. The 
follow up measurements reveal behaviour consistent with dust as a 
scatterer, i.e. an increase as /!“' between g' and F before a less steep 
increase in polarisation for the 425SP filter; this is shown in Figure 
10. The polarisation angle in all three bands is the same within the 
error of the measurements. 
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Clarke & Schwarz (1984); Schwarz (1986) were able to show 
that circumstellar dust is not the only mechanism reponsible for 
polarisation of red giants. Harrington (1969), in reference to Mira 
variables, was the first to propose that Rayleigh scattering from 
the photosphere could be responsible for the observed polarisa¬ 
tion. This mechanism was initially ruled out by Dyck & Jennings 
(1971) on the basis that the wavelength dependence was incon¬ 
sistent with the dependence of Rayleight scattering. It wasn’t 
until spectropolarimetry of supergiant Betelgeuse (a Ori) revealed 
polarisation that changed across, and was correlated with, spec¬ 
tral lines of, in particular, TiO (Clarke & Schwarz 1984; Schwarz 
1986), and that the relationship seen was an artifact of broad¬ 
band polarimetry that this mechanism was accepted. This mecha¬ 
nism of Rayleigh or Thompson scattering from hot-spots or con¬ 
vection cells in the photosphere is able to explain polarisation that 
varies regularly (in regular variables) on time scales of months to 
years, whereas alternate mechanisms such as non-radial pulsations, 
or an equator to pole temperature gradient cannot (Schwarz 1986). 
If there are multiple hot-spots a polarisation angle that varies fairly 
smoothly across bands can also be expected, however that is not 
the case here. Hot-spots as the origin of polarisation doesn’t pre¬ 
clude the presence of additional scattering regions as well. In fact 
changes in polarisation across TiO lines point to TiO as an absorber 
or scatterer in the outer photosphere of highly polarised stars. It 
seems reasonable to expect that should hot-spots be the cause of 
polarisation in 6 Oph, then molecules like water, be they present 
in an extended cool photosphere or in a MOLsphere would also 
have spectral lines with a polarisation different to the surrounding 
spectral continuum. 

Starspots are also a possibility, with Saar & Huovelin (1993) 
having calculated maximum polarisations for late giants from this 
mechanism of a few hundred ppm. Being evolved stars red giants 
are generally thought of as being inactive. However, a wavelet anal¬ 
ysis carried out by Hedges et al. (2013) found evidence of starspots 
in 14 out of a subsample of 416 red giants observed by the Kepler 
satellite. Very recently significant magnetic activity has been de¬ 
tected in Mira (Vlemmings et al. 2015), and the RS CVn K giant 
V* XX Tri is famous for being observed with a ‘superspot’ cov¬ 
ering 11% of its surface area (Kunstler et al. 2015). Hedges et al. 
(2013) analysis revealed features that could persist for a few hun¬ 
dred days. Even so, this is not the norm, and seems the most un¬ 
likely candidate for the polarisation seen in 6 Oph. The consistency 
seen in our two g' measurements a year apart would put 6 Oph at the 
extreme end of the range with regard to the persistence of such fea¬ 
tures. The maximum polarisations calculated by Saar & Huovelin 
(1993) rely on a filling factor of 24%, and something approaching 
this level would have to be present to see a ~500 ppm polarisation. 
This doesn’t fit for a star not known to be magnetically active. 

Furthermore, similar degrees of polarisation are calculated for 
g' by Leroy (1990) due to differential saturation of spectral lines 
in magnetic fields. Leroy (1990) was most interested in K dwarfs, 
but also looked at a K giant. That analysis shows that polarisation 
due to differential saturation should be greatest in the 425SP band 
- which is not the case for S Oph - a worked example spectrum for 
Arcturus (BS 5340, K1.5III) with a hypothetical 500 Gauss field 
is shown. Arcturus was observed as part of the PlanetPol survey, 
recording 6.3±1.6 ppm; whereas in PlanetPol’s waveband the cal¬ 
culation yields around 150 ppm. If this mechanism were primarily 
responsible for the polarisation we see in M and K giants we might 
expect the ratio of median polarisations determined by HIPPI and 
PlanetPol to be similar for M and K spechal types, which is not 
what we see in Table 6. 


5 CONCLUSIONS 

Fifty of the brightest stars in the southern hemisphere have been 
observed for their linear polarisation. We have investigated these 
alongside 49 northern stars from the PlanetPol survey. These two 
surveys are far more sensitive than previous stellar polarimetric sur¬ 
veys. Prior to this work, aside from the four classical Be stars, none 
of the stars had identified intrinsic polarisation assigned to them. 
We have now attributed polarigenic mechanisms to around half of 
the stars in the survey. In particular B type stars have been identi¬ 
fied as being polarised through a variety of mechanisms, whilst late 
giants also show a greater propensity to be intrinsically polarised 
than demonstrated in previous studies, mostly as a consequence of 
circumstellar dust. Intrinsic polarisation was also identified in an 
Ap star and in debris disk systems. 

Amongst the intrinsically polarised B stars was the close bi¬ 
nary system Spica, which demonstrated variability from repeat 
measurements. Of the four Be stars observed, three have been re¬ 
solved by interferometry. For each of these we measured a polarisa¬ 
tion angle aligned perpendicular to the rotational axis. We recorded 
the highest degree of polarisation on record for the Be star Achemar 
in September of 2014; this indicates it was particularly active at 
this time. It has since declined in polarisation. Other B type stars 
were also found to possess high levels of polarisation. Since most 
of the B stars rotate at well below break-up velocity (Regulus in 
the PlanetPol survey being an exception) this cannot be a result of 
rotational oblateness. Instead we advance the hypothesis that what 
we are observing is low level Be-like behaviour in normal B stars. 

There are only a small number of debris disk systems in the 
combined catalogue of the two surveys, but these show a slight ten¬ 
dency to be more polarised than other nearby stars. The circumbi- 
nary debris disk system e Sgr appears to have a polarisation signa¬ 
ture consistent with a companion-induced asymmetry in light scat¬ 
tered from the disk. 

Many of the M and K giants in the surveys are clearly intrinsi¬ 
cally polarised with cooler and more luminous stars having higher 
polarisations. Infrared excess is not correlated with the degree of 
polarisation in our sample. However circumstellar dust distributed 
by an asymmetric stellar wind may be responsible for the trends 
we see. As might an increased scattering density from dust closer 
to the photosphere in cooler stars. However the dust free M giant S 
Oph also displays a variable and high degree of polarisation. This 
behaviour may be associated with Rayleigh scattering from stellar 
hot-spots, or less likely star spots. 

The plethora of polarised stars identified has profound conse¬ 
quences for polarimetric studies of the local interstellar medium, 
and for searching out other polarimetric effects nearby. It is com¬ 
mon practice to determine a value for the interstellar polarisation of 
an object thought to have some interesting intrinsic effect by tak¬ 
ing an average of the measured polarisations of nearby stars. Stars 
are selected for this purpose on the basis that they are likely to be 
unpolarised. Our study shows that if one seeks a value for the in¬ 
terstellar polarisation at the ppm level then, in particular, no B type 
stars nor M or K giants should be selected. Ap stars and debris disk 
systems should also be avoided. As these are some of the brightest 
stars at any given distance and otherwise likely to be selected this 
is an important finding. 

When we discount the intrinsically polarised stars from our 
sample and that of the PlanetPol survey, we find that interstellar 
polarisation is stronger in southern stars. Away from the Galactic 
equator in the northern sky stars are polarised at ~2 xlO^’ pc“*. 
At distances between 10 and 30 pc in the southern sky the rate 
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of polarisation can be as much as ~2 xlO“^ pc“*. Beyond 30 pc 
there appears to be less of a polarising effect from the interstellar 
medium. The polarisation angles of stars just to the Galactic north 
are aligned to the Galactic equator. A dustier interstellar medium 
to the south (in the direction of the centre of the galactic plane) 
appears primarily responsible for the higher degrees of interstellar 
polarisation observed in southern stars within 100 pc. 

The very high precision of HIPPI in fractional polarisation has 
allowed us to uncover the polarimetric behaviour of some of the 
best known stars in the night sky. We have demonstrated HIPPFs 
potential in many areas of stellar polarimetry, producing more 
precise measurements than have previously been possible and in¬ 
deed uncovering the polarigenic behaviour of some stars previously 
searched for but not found. HIPPPs simplicity and economical fab¬ 
rication belies its great value. 
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